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SAMPLED-DATA SYSTEMS WITH EXTRAPOLATING DEVICES* 


la. Z. Tsypkin 
(Moscow) 


Automatic sampled-data systems containing extrapolating 
devices are investigated. Equations for such systems, describing 
the processes at any moment of time, are given. The method of 
analyzing such systems is illustrated by examples, 


INTRODUCTION 


Recent years have been characterized by the widespread use of discrete methods of measuring,. transmitting 
and transforming physical quantities. Discrete methods of measurement give the possibility of measuring physi- 
cal quantities without using devices the states of which are characterized by these quantities, and also allow the 

measurement of those quantfties whose definition involves the passage of time. 


Fig. 1. 


Discrete methods of transmitting and transforming are distinguished by their high accuracy and impervious- 
ness to noise, and allow one and the same channel to be used for the transmission of a large number of data. These 
methods are used in automatic sampled-data systems and in automatic systems with digital computers. 


In many cases, the use of such systems gives rise to the necessity of reproducing a continuous function from 
data transmitted only at equally-spaced intervals of time, Tp This problem leads to the reconstruction, in the 
time interval nT), < t < (n+ 1)Tp, of the continuous function from its values at the previous moments of time, 
t= mT), wherem=n~1,n—2,.. .»n~k, that is, to the problem of extrapolation, 


Since discrete data gmerally contain random components, it is necessary, in addition to extrapolation, to 
smooth or adjust the discrete data. 


The investigation of sampled-data systems, designed to transform discrete data to continuous data,was expounded 
in works [1-3]. In this paper, we consider automatic sampled-data systems in which an extrapolating device (ED), 


* This paper contains one section of the talk, “Certain questions in the theory of discrete systems,” presented at 
the Symposium on the Applications of Computers in Automatic Systems, held at Atlantic City in October, 1957. 


i. €., a transformer of discrete data to continuous data, ts one of the system elements. A block diagram of such 
a system fs given in Fig. 1. The figure also gives the character of the variation of the output of each system ele- 
ment, The darkened segments show the computation of the output z of the continuous portion of the system 
from the input quantity f. 


The output of the data-sampling (DS) element (discrete data) is fed to the input of ED, which transforms 
it into a continuous quantity which changes, between sampling times, according to a definite law. Generally, 
this law is taken as a polynomial of degree s. Such an ED will henceforth be called an s~order ED, or transformer. 
A zero-order corresponds to a fixed element, which retains a discrete value until the appearance of the next 
discrete value. A first-order ED is implemented by a linear extrapolator, etc. The output of the ED acts on the 
continuous portion of the system (CP). 


Since the extrapolation under consideration is intimately related to the known character of the processes 
during the repitition periods, 1. e., between the discrete sampling times, the equations for a system with an extra~ 
polating device must describe the processes at any moment of time. 


1. Structure and Equations of the Extrapolating Device 


We begin by clarifying the structure of an ED. 


If the law for transforming the discrete data to continuous data is restricted to be a polynomial of degree 
s, then the output of an s-order ED may be represented by the equation 


8 
y(n, e] = a, (1) 


v=0 


where 0 < € <1 Is the relative time measured from the beginning of each interval t = n, and the coefficients 
a, {0} change from interval to interval. In the general case, the coefficients satisfy the difference equation 


[nm] = po 0] —(1 — po) — pI. (2) 
TABLE 1 
Extrapolation law 
[n —v, 0 
0] 
2 Newton's Formula a, = > (— 
pev 
— Sterling number 
k 
Method of least = 4 've 0}, 
3 
linear extrapolation is si —v, 0 
( po ) ay [n} 2“ 2u+ 1) O] 


Here pp is a constant multiplier, ay ), gi”) are constants, y,[n, 0] is the input and a, [n] is the output. 
For pp = 0 or pp = 1, the magnitude of a, [n] depends only on, respectively, the output or input quantity. For 
0 < py < 1, the magnitude of a.,,[n] depends both on the input and the output, and the degree of this dependency 
is determined by the magnitude of pp. 


ex, ot,/nje* ex,/nje" 


Fig. 2. 


The law of variation of the coefficients a,,[n] depends on the choice of the extrapolation law. Different 
methods of extrapolation are quite well described in the specialized mathematical literature (cf., for example, 
[4}). Certain particular cases of laws of variation of the coefficients «,,[n] are given in Table 1. 


From Formulas (1) and (2), it is possible to deter- 
mine the structure of the ED. 


Each term a,[n] in (1), during the interval 
when € varies from 0 to 1, determines a “graduated 
pulse” (Fig. 2), 1. e., a pulse of parabolic shape, with 
height equal to a,[n], which varies from interval to 
interval according to the law defined by difference 
equation (2). 

The graduated pulses can be formed by supplying 
instantaneous “square” pulses, a,,[n}, to anelement 
whose transfer function equals the representation of the 
graduated pulse [5). 


Thus, the structure of an ED will have the form shown in Fig. 3. 
The discrete data are supplied to the input of discrete filters 1, described by Eq. (2). The filter output, 


also discrete, is fed to the input of forming elements 2. The graduated pulses at the output of these elements, 
when added, form, in each interval, the continuous curve y [n,€]= y(t). 


To determine the transfer function of the extrapolator, we subject difference equation (2) to the discrete 
Laplace transform [6, 7]: 


(q, = In, e] = D{/[n, e}}, 
n=0 


where € = 0. 


Using the notation A, (q) = D {a,(n]}, 0) = D{y{n, 0}} , we obtain the equation of the dis- 
crete filter with respect to the representation: 


A, (q) = Ka, (q) 0), (3) 
where 
Kay (q) = (4) 


1+ (1 — po) 
p= 


is the transfer function of the discrete filter which, for pp < 1, may be considered as a sampled-data system with 
feedback. 
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Fig. 6. 
A discrete filter may be realized in the form of a circuit with lag elements and amplifiers (Fig. 4a). A 
more detailed layout of the circuit is given in Fig. 4b [6). 


It is also possible to use, as a discrete filter, a digital computer (DC), the program for which is determined 
by difference equation (2) or by transfer function (4) [6]. 


By definition, the transfer function of the forming elements must coincide with the transform of the gradua~ 
ted pulses, i. e., 


T v 
v) v—p) 
Kp (p) = dt = >» (5) 
where 


Ti} 
a b 
Fig. 4. 
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is a factorial polynomial. 


It follows from (5) that the forming elements must contain lag elements, amplifiers and integrators. The 
schematic for one of the forming elements, for example, the v~th, is given in Fig. 5. With increasing v , the 
number of integrators increases. In practice, it is not necessary to use different integrators for each v-th (v = 
= 0,1, 2,..., 8) forming element, It is possible to use the integrators of the s-th forming element for all form- 
ing elements, as is shown in Fig. 6. Here, the sign £ denotes a summing element. 


Thus, an ED consists of two parts, a discrete part (discrete filters) and a continuous part (integrators) [2]. 


2. Equation of a Closed Sampled~-Data System with an Extrapolating Device 


A schematic for a closed sampled-data system with an ED may take the form given in Fig. 7. In the total 
continuous portion of the system (total CP) are now included in the integrators of the forming elements (FE). If 
we denote by W_-(p) the transfer function of the continuous portion of the system, then the transfer function of 
each channel, consisting of the continuous portion and of a forming element, will equal 


W,(p) = Key (p) (Pp) 


or, in dimensionless quantities, 


Fig. 7. 


This relationship follows from the similarity theorem, according to which if Wy(p) = L{w,(t)} , then 


We substitute into (7) the value for Kzy(p), first replacing p in it by + , We get the result 
P 


v—p) 
pmo? 


i 
Here, we took W,.(q) = ) 
Tp 
For determining Wy, (q), it {s convenient to use Tables 2 and 3, which establish the relationships between 
W(q) and W*(q). The derivations of these relationships may be found in, for example, [7, 9, 11]. 


We shall assume initially that Wc (q) is given in analytic form. On the basis of Formula(3)in Table 2, we 
can also find W- (q, €) in analytic form. With the use of Expressions 1 and 5 of Table 3, we can write the trans- 


fer function W,(q , €) corresponding to (8), in the form 


| 
1 q q 
Wy = K Ww. (7) 
(a) = Me (zt) 
where q = pT). 
Oy | total CP 
a 
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TABLE 2 


we w@ Remarks 
w |n, = ¢] = w (n + €) 
= Kuss: i 
ttern 
w a for t=n+e 
W* (q, = \ 
2 |= > efatenine W (q + Qnir) 0 Transfer function 
r= —o of the continuous 
portion 
W (q)= Poles W (a). 
Cc =0 ipl 
W*(q,¢) = > > x s ae 
3 v=0 p=0 => > ple r, (v= i. 
\ 
Here 
1 ry 
TABLE 3 
(q) (q, e) 
k — integer 
(y <1) w*(q,e—y) (rset) 
©) 
3 qW (q) 
lim q’W (q) =0 ae” 
4 € 1 
5 \w (q, ede += —\ Wa, eae 
€ 4 1 1 
n=0 0 
i v times 
—W 
6 (q) where 
i 
0 
and . 
(9, ©) = W"(q, ©) 


” 
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Wy (q, = e) dev) — 


owe 


(Wola, e) 4 
0 


v—p+1 times 


¥ 1 


(v— »)! 


ef— 1 


Win-1) (q, de — 


0 


where 


i 


Win (a, = \ Wena (a, de + 


1 
(q, de 
Woo) (q, 8) 


The computation of the integrals with respect to €, which enter into (9), presents no difficulty, since € 
enters either as a multiplier or as an exponent, 


The transfer function of each channel, taking the discrete filter (DF) into account, is found by the simple 
cross-multiplication of the transfer functions (4) and (9), 1. e., K%,,(q) Wy (q, €). 


We can now easily determine the transfer function of the open-loop sampled-data system with an ED; 
W" = Kav (a) Ws (4, (10) 
v=0 


The equation for the closed-loop sampled-data system with an ED is written according to a formula known 
from the theory of discontinuous regulation 
— 
= 
in the form 
8 


(q, =—— F* (q, (11) 
i+ >) 
v=0 


: Kay (2) Wy (9 ©) 
K,(q, 8) =—>, (12) 
1+ > ©) 


is the transfer function of the closed-loop sampled-data system with an ED. If we replace q in(12) by jw = 
= jw Tp, we obtain the frequency characteristic of the closed-loop sampled-data system. 


Particular cases subsumed under Eq. (12) are the equation for a sampled-data system with an indexer (for 
s = 0) (3, 5}, and the equation of a system for transforming discrete data into continuous data (for W-(q) = 1) [2]. 
3. Different Forms of Presenting the Transfer Functions of Open-Looped Systems 


In those cases when the transfer function of the invariable continuous portion W_(q) {s unknown or com- 
plicated, it is possible to use experimentally obtained impulse or frequency characteristics to determine the trans~ 
fer function W*(q, €), or the frequency characteristic W* (jw ,€). 


Let the impulse characteristic of the continuous 


wf ty) w(t) bu/n9) portion w(t) have the form represented in Fig. 8a. 
As follows from (8), the impulse characteristic 
A ; for each channel, exclusive of the discrete filters, is 
” 
a b 


written in the form 


Fig. 8.. 


t 
w, = TY dit? — (3) 


Integrating graphically v times(v = 1, 2,..., 8) the impulse characteristic , and making use of Formula(13), 
we find w,,(t) (Fig. 8b). Since w,,(t) grows with increasing Tf, the following formula should be used to determine 
Wy (a. €) 


ef at w, [0, 
4 ’ 


Ws (q, = A* w, [n, + (14) 


v=1 


which follows from Formula(1)of Table 2 if the theorem on the transform of differences [7] is used. There may 
always be found a k such that 


lim A* w, [n, e] =0, (15) 


/ 


and in practice, (14) may be limited to a finite number of terms,* The expression thus found for Wy (q, €) are 
substituted into Eq. (11). 


* A similar formula, for k = 1 and q = fw, was used in [8]. 


| 


In the case where it {s difficult to find an analytic expression for Wi,(q, €), given the frequency characteristic 
of the continuous portion Wc(jw), it is easy to find the frequency characteristics Wy(jw,¢), and thereafter to 
use the frequency method [9] to compute the sampled-data system, 


According to Formula(2)of Table 2 we have, forq=j@, 


Wy (7G, 2) = (G + = 


= + 2ar)] We (@ + (16) 
where Kp, = (ie 
. Substituting Wi Go, €) into (10) for q = jw, and changing the order of summation, we obtain, taking perio- 
dicity into account, 
K  (J@) = (@ + 2nr)), (17) 


with the expression for W°(jw, €) in the form 


2) = (@ + x 
v=0 


As a consequency of Formula (18) we obtain the following method for constructing the frequency charac- 
tersitic of an open-looped sampled-data system with an ED. 


Changing the numerical frequency characteristic W(jw) by the relationship w = wT 
1 its modulus by 1/Tp, we obtain Wc. (j w). Rotating the radius-vector of Wc (jw) 
(13), angle @€ (€ = const), we find 
ei W (ja). (19) 


Then we construct the curve 
>) Kpy (78). (20) 
v=0 


Cross-multiplying curves (19) and (20), we find 
v=0 


Using well-known methods [9], we obtain W°(jw, €) from this curye in accordance with Formula (18). 


4. Determination of the Processes in Closed-Loop Systems with Extrapolatin 
vices 


Methods known from the theory of discontinuous regulation may be applied to Eq. (12) for the investiga- 
tion of sampled-data systems, 
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If W* (q, €) is found in analytic form, K* (qd, €) may be given in the form 


ca 
where H*(q, €) and G *(q) are polynomials in e4. 
The impulse characteristic of the closed-loop system, k, [n, €], is found from the decomposition formula 
[7] 


n, e] = >, ev 23 
ke[n, tay (23) 
Here Gy are the poles of K‘, (q, €) which, for simplicity, we shall assume are simple poles. 
If W* (jw, €) is found graphically then, by means of the method described in [9], it is possible to find the 
frequency characteristic of the closed-loop system, or its real part Re Ke(j@,¢€). Then, from the formula 


k In, =2 (Re e) cos @nda@ (24) 


it is possible to determine the impulse characteristic of the closed-loop sampled-data system. 
The process for an arbitrary impulse f{[n, €] is found from the superposition formula 


z[n,e)= >) k jm, e}f[n—m, 0}. (25) 


5. ‘Example of the Computation of Systems with Extrapolating Devices 


As an example, we consider an automatic system with a linear ED, 1. e., with a first-order transformer. 
The block schematic of this system is shown in Fig. 1.* We shall assume that only the frequency characteristic 
of the continuous portion of the system without the extrapolator is known (Fig. 9). 


-1 - — D 
a8 15 157 
IN 4-1 
fas +2 
+5 
4-7 
4-8 “7 
9 
4-19 


Fig. 9. Fig. 10. 
* A similar system, but with much simpler assumptions, was studied experimentally on a model in [10]. 


n 
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The extrapolation law is determined by Eq. (1) for s = 1, where 


aan} (Sts Im. 0} + [n — 1.0] — ys — 2.0), 
[n] = (ys [n, 0] — ys [n —2.0}). 


These values are found by the method of least squares from the condition that the best straight line be 
passed through the previous three points (Formula(3)of Table 1 for k = 3). 


According to Formula (4), the frequency characteristics of the discrete filters will be 


(26) 
(i@) = (1 — 


The frequency characteristics of the forming links, 
according to Formula (5), with p = £2 » will have the 
form: 


tine) = 7, 


2 


/ On the plane of the frequency characteristic 
‘ Wc (f@) we construct the curve from the expression 
Fig. 12 by Wc (fm) and extending this periodically, we find 
Roe W*(j@, 0). Figure 10 shows the curve for W*(jw, 0) 
for T), = 2 sec and ky = 1, where kg fs the total gain; the dotted curve gives the frequency characteristic for the 
open-looped system with a zero-order extrapolator (a,[n] = 0). 


As is clear from a comparison of these characteristics, the use of a linear extrapolator lowers the stability 
of the system. From W*(jw, 0), using ordinary graphic methods [9], the real part of the frequency characteristic 
of the closed- system ReK¢(jw, 0) is constructed. This is shown in Fig. 11, where the open-loop gain is 
taken to be ky = °/s. By approximating Re K¢(jW, 0) by the sum of elementary trapezoids, we find the impulse 
characteristic of the closed-loop system k,([n, €] as the sum of the impulse characteristics corresponding to the 
individual trapezoids (Fig. 12a). The process for jump impulses is determined according to Formula (24) using 
simple summation, and has the form shown in Fig. 12b. 
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A WAY OF FORMING TRANSFER FUNCTIONS OF SAMPLED-DATA CONTROL 
SYSTEMS WITH EXTRAPOLATING DEVICES 


S. S. Ermakov and E. M. Esipovich 


(Moscow) 


The paper deals with a way of forming condittonal transfer 
functions of extrapolating devices for digital-to-analog conver- 
sion, The transfer function expression depends on the shape of 
the input pulses. 


1. Characteristics of Considered Control System 


There exist control systems which contain a section (pulse unit) which converts a continuous function of 
time into a series of its values equally spaced in time, and which are characterized by the fact that the object 
being controlled is under their continuous influence. This last condition distinguishes these systems from pure 
pulse systems in which the object being controlled is acted upon periodically, during certain intervals of time, 
beginning with time nT», where n= 0,1, 2,.. ., and during the rest of the time is not subjected to any control- 
ling action. 


One way of achieving continuous control over an object, with systems containing the above-mentioned 
pulse unit, having a small build-up,is by digital-to-analog (or sectionally~continuous) conversion using extra- 
polating devices (extrapolators). 


The present work sets forth a method of forming 
the transfer functions of extrapolating devices, which 
permits the application of the existing theory of sam~ 
pled-data control to the analysis and synthesis of sys~ 
tems containing these devices. 


Figure 1 shows a block schematic of the type of 
control system being considered, where PU is the pulse 
unit, ED is the extrapolating device, LS is the linear 
(continuous) section of the system, including the object 
being controlled. 


Discrete values Au (nTp] of the continuous function 


Au (t) = (t) — (1) 


permit the construction of the sectionally-continuous (or continuous) function Aug(t), similar in a specific sense 
to the function Au(t), which we are considering, and thus fill in the gap in our knowledge of the two functions 
being compared, ujn(t) and Ugy,(t). This operation is carried out by the extrapolating device ED, which thus 
predicts the nature of the change in the function Au(t) between the times nT, when the measurement {s made. 


Let us state one of the possible methods of extrapolation. 


q 
| 
Fig. 1. 
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Tp<t<(n+1)T, (2) 


let us draw sections of straight lines (Fig. 2) in such a way that the sum of the squares of the inclinations of each 
straight line from the k of the preceding points {s a minimum. 


elt) 


8 


4. 


(2-2), 


Fig. 2. 


The equations of the indicated straight lines can be easily derived in a general form, but we shall limit 
ourselves to the simplest example of linear extrapolation from three preceding measurements. The method pre- 
sented below can be applied to more complicated cases. 


Direct computations, which result in the finding of the minimum of a function with several variables, lead, 
‘in the case of extrapolation from three preceding measurements, to the following expression: 


= — Ian) + + — (3) 


Tn / p], [(n 1) Ty), / [(n — 2) T>) (4) 


are the measured values of the function f(t) = Au(t), from which the segment of straight line f,(t) which pre- 
dicts the behavior of the function f(t) on Interval (2) is constructed. 


Let us call the magnitude f, [nT], corresponding to the beginning of Interval (2) (Fig. 2), the average 
value of the function, and the ne i fel(n+ 1)T, ], corresponding to the end of Interval (2), the predicted 
value, The introduced notation permits one to distinguish between the beginning and end of each piece of the 
sectionally~continuous function (3). The discontinuities in the extrapolated function f,(t), constructed in accord- 
ance with (3), have an undesirable effect on the dynamics of the control system. These discontinuities can be 
avoided by changing the method of extrapolation. 


One of the possible methods of obtaining a continuous extrapolated function f,(t) is shown in Fig. 3. Here 
the position of the predicted point B is determined from the three preceding measurements; f [(n— 1)Tp], 
f((n— 2)T,) and f ((n- 3)Tp}]. At time nTp the value f[nTp] is fed to the extrapolator. The presence e this 
value together with the values f [(n— 1)Tp jand f 2)Tp] permits, at this same time nT», the 
determination of the position of the following predicted point D. It is therefore possible to move in the direc- 
tion BD while avoiding the jump from section AB to section CD. The equation of the straight-line section BD 
is determined, in our example, from the four preceding values of the input variable f(t). 


>? In the intervals 
Not 
=~ 
In Eq. (3) 


| 


8 


1 


Fig. 3. 
It follows from what has been said that the graph of the output variable from the extrapolator with the con- 


tinuous output (Fig. 3) can be obtained from the consecutive connection of the predicted points of all the inter- 
vals. The corresponding analytical expressions are tabulated in the Appendix. 


If the extrapolating device is looked upon as a section of the control system, then the extrapolating equa- 
tions tabulated in the Appendix (among them Eq. (3)), establish the dependence of the output variable f,(t) on 
the input f(t). This dependence is expressed by means of a differential equation. 


The other characteristic of the extrapolator is the fact that the output f,(t) is determined only by the values 


re= flv Tp] = 0, 1, 2,.. .) and does not depend on the form of the function f(t) between the points of measure- 
ment, yT,. It is possible to construct many functions f(t), coinciding at the times » T,which correspond to 
lead one and the same function f,(t) at the output of the extrapolator. It is, therefore, not possible to form the trans- 
‘ fer function of the extrapolating device (as a relationship between the input and output) using methods based on 
the theory of continuous control. 
) 2. Derivation of the Transfer Function of the Extrapolator 
The extrapolator can be looked upon as consisting of the pulse unit PU and the linear part with the trans- 
fer function W,(p) (Fig. 4). The pulse unit converts a continuous function f(t) into a multistep one f[v Tp). 
In accordance with the theory of sampled-data control 
) [1], such a unit is characterized by an amplification 
ac il fh. factor kj = 1 build-up y = 1. 
pre- | PU tft) 
We shall consider the multistep function f[v Tp] 
Fig. 4 to be the input of the section, representing the linear 
ig. 4. portion of the extrapolator, and the function f(t) to ; 
ted be the output of this section. Equation (3) expresses the relationship between them. 
Ss Taking the Laplace transform of both sides of Eq. (3): 
be 
t—aT 
fet) = | Un — dt + 
this P 
+\ (fm + fama) (5) 


f(t), ft) | 

| 

1 
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In view of the fact that the integrand on the right side of (5) changes from interval to interval, we shall 
_ break up the integrals in (5) into a sum of integrals: 


(n+1)Tp 


(t)e-*t dt = >} 


P (fn — fn—2) ™ dt + 


P 
(n+1)Tp 


+¥ + 2fn—1— dt. (6) 


Utilizing the fact that the function f, = f[v Tp] is constant over the limits of every integral in (2), we 
transform (6) to: 


(n+1)Tp 
felt)e- dt = — \ (t — nT) dt + 
0 Pp aT 


(m+1) Tp 


Simple calculations, with due attention to the designations in (4) yield: 


\ fe (t) dt = 5 pr, —e pT ye x 
0 


x >. {f — f — 2) Ty]} + x 


{5f + 2f [(n —1) Tp] — f — 2) 


Assuming, as usual, that f(t) = 0 for t< 0, it can be shown that 


> — a) Tp} = f [nT 


Equation (9) permits the simplification of Expression (8): 


i —pT. —2pT. i—e —pT 


| 
(9) 
n=0 n=0 
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On the other hand, we find what must be the expression for the input variable f (nT p): 


0 n=0 
—pnT 


P n=0 


Comparing (10) and (11), we find the transfer function of the linear part of the extrapolator: 


We(p) = opr, (! — — pT (4 + P79) + 


(5 + ny, 


as a relationship between the expressions for the output and input values, 
Let us call Operator (12) the conditional transfer function. Further investigations of the system are carried 
out following well-known rules of sampled-data control. 


The question of the formation of transfer functions of extrapolating devices has been also considered using 


a different approach, by Ia. Z. Tsypkin [3]. The author represents the extrapolator by acomplex pulse unit which 
forms pulses of a definite shape (see [2]. 


It can be shown that the pulse transfer function [4] 


We (q+ 2nmi) + (13) 


(where q = pTp), built up by the method of conditional transfer functions, is identical to the corresponding trans~ 
fer function of the extrapolator derived by la. Z. Tsypkin [3]. 


For this we must make use of Relationship (8), given in [2], and take into account the periodicity of the 
transfer functions of discrete filters introduced, in accordance with [3], into the circuit of the extrapolator. Then 
it is not difficult to transform the transfer function, derived by Ia. Z. Tsypkin, into the form of (13). It is also 
possible to achieve the reverse transformation into a conditional transfer function. 


Calculations have shown that the use of the indicated transfer function is convenient in those cases where 
it is necessary to examine the effect of different methods of extrapolation on the operation of control systems. 
The application of conditional transfer functions shortens calculations in cases where one {fs interested in the 


choice of the optimum parameters of the linear part of the system without changing, at the same time, the con- 
stitution of the extrapolator. 


SUMMARY 


In the study of the dynamics of a control system containing a series-connected pulse unit, having infinitely 
small buildup (y —> 0), and an extrapolating device, these sections can be replaced by a pulse unit forming rec- 
tangular pulses having buildup y = ba and the linear (continuous) part of the extrapolating device, also connected 
in series. 


The transfer function of the linear part of the extrapolating device can be found by means of the usual (not 
discrete) Laplace transforms. This transfer function expresses the connection between the expressions for the input 
and output functions only for a particular form of input — a series of rectangular pulses having buildup y = 1. 


Hence such a function may be called a conditions] transfer function. 


The above material allows us to assume that the method of conditional transfer functions (in the indicated 
sense) 1s applicable when the input consists of a series of pulses of arbitrary, known form. 


APPENDIX 


The table lists formulas relating the output variable f,(T) of the linear extrapolating device to the discrete 
values of the input variable f(T), as well as the corresponding conditional transfer functions W¢ (q). 


The time T = t/Tp is expressed in terms of the interval time Tp. Hence, in particular, q = PTp- 


Method of 
extrapola- Extrapolation formula for Interval (2) Conditional transfer 
(n = 0, 1, +) function We(q) 
tion Ss 
fe ®) = 
i 
ES-2 fe(®) = {f — — 2) +f Hi 
i 
— ge? 
+ 2/ — 1] —f In —2}} 
4 
ES-4 @ = +f 1] — Tog (B+ — x 
— f 2] —3f [n—3]} + x(1 — — ge%) + 
— 267) 
i 
ES-5 = + — 1] — (2 + + + 
—f |n— 3] — 2f [n—4]} + — ge) + 
4 
i 
+ 2f |n—1] —f [n—2 
4 
EC-4 felt) {44 fn] — 3/ — gq (4 — 22%) x 
— 3f [n — 2] + 2f [In —3]}} + x(1— — +. 


Method of Extrapolation formula for Interval (2) Conditional transfer 
extrapolation (n =0, 1,.2,...) function W, (q) 


EC-5 


sts i 
felt) = {2f [n] —f — (2 +e — x 


— 3f 3] +4 + — 
+ + 5/ In — 2} 4.2 — 


The following designations are introduced in the table; ES~-3 is an extrapolator with a sectionally-contin- 
uous output function f.(0), determined from the three preceding values of f[v](v =n, 2-1, n— 2) of the input 
function f(t) (see Fig. 2); EC-4 is an extrapolator with a continuous output function f,(t), determined from the 
four preceding values f(v](v =n, n—1, n— 2, n— 8) of the input function f(T) (see Fig. 3), etc. 

LITERATURE CITED 


[1] Ia. Z. Tsypkin, Transient and Steady-State Processes in Pulse Circuits [in Russian] (Gosenergoizdat, 
1951). 


[2] Ia. Z. Tsypkin, “Concerning the calculation of the pulse-shapes in intermittent-control systerns,* 
Automation and Remote Control (USSR) 16, 5 (1955). 


(3) Ia. Z. Tsypkin, "Sampled-data systems with extrapolating devices,” Automation and Remote Con- 
trol (USSR) 19, 5 (1958).* 


[4] Ia. Z. Tsypkin, “A frequency method of analysis of intermittent-control systems,” Automation and 
Remote Control (USSR) 14, 1 (1953). 


Received November 1, 1957 


*p. 383, this translation, 


AN ANALYTICAL METHOD OF SYNTHESIS OF LINEAR CONTROL SYSTEMS 
IN THE PRESENCE OF INTERFERENCE AND WITH A SPECIFIED DYNAMIC ACCURACY 


K. I. Kurakin 


(Moscow) 


An analytical method of approximation of transcendental 
transfer functions of automatic control systems is proposed. The 
functions were obtained with the help of fraction-rational func- 
tions [10, 14]. The practical applications of the method are {l- 
lustrated by several examples. 


INTRODUCTION 


The raising of the requirements of modern technology for speed, dynamic accuracy, and quality of auto- 
matic control systems, providing separation of signal from noise, brings out the necessity for the development 
of an engineering method for the calculation of these systems, with a view to obtaining optimum designs and an 
accurate evaluation of their working characteristics, This is particularly important in the fields of radar track- 
ing systems and automatic control systems. In a general case there appears at the input of the indicated systems 
a useful signal 


Bi (t) = + ate (t), ay 
k=l 


composed of a random signal s(t) and a series of known functions of time f;,(t),.. ., fy(t), as well as a random 
noise signal By (t) = n(t). Moreover, s(t) and n(t) may be stationary or nonstationary random signals, while the 
coefficients ay, . . ., an may be arbitrary, unknown quantities or random variables with a given mutual correla~ 
tion 


E [aaj] = yij, bi (2) 


The theory of the interpolation and extrapolation of stationary random sequences was originally developed 
by A. N. Kolmogorov [1]. Subsequently, the problem of the optimization of a linear, dynamic system for the 
case when the input consists solely of continuous stationary random functions s(t) and n(t), was solved in [2). 


It was shown in [3-6] that for automatic control systems it is often more expedient to examine and establish 
the statistical properties of the derivatives of the useful signal s(t). The resulting optimum solutions, obtained 
in this case, are related to a class of linear, astatic systems which have practical applications. The statistical 
theory of optimization, developed by Wiener, was generalized in [7] for the case when the signal and noise are 
both nonstationary random functions, The optimum pulse transfer function of the system w(t,r), for a nonsta~ 
tionary input, represents a linear system with variable parameters, its response being determined by the super- 
position integral 
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Bo (t) = =) (2) + Bn de. (3) 


A method of constructing an optimum system, with constant (in time) parameters, which ensured, in a 
specific sense, the filtering of nonstationary random functions, was proposed in [8, 9]. The theory of the opti- 
mization of linear dynamic systems was further advanced in paper [10], for the case when the observation time 
T is finite, s(t) and n(t) are the stationary, random functions, and the mathematical expectation of the function 
s(t) is unknown, but is determined on the interval T by a function of time in the form of a polynomial of degree 
r. Davis [11] generalized the problem, formulated in [10], for the case of nonstationary random functions s(t) 
and n(t). V. S. Pugachev [12] derived the general, necessary and sufficient condition for a minimum in the 
mean square error of a dynamic system, this condition being the source of all the known equations of optimum, 
dynamic systems studied in the above categories. Among the shortcomings of papers [10, 11] fs, firstly, the use- 
lessness of the final equations, from the point of view of their physical realizability [13] for practical calcula~ 
tions, An attempt was made in [14] to develop a method for determining the optimum pulse transfer function, 
stated in [10], for its application in the calculation of linear tracking systems. By means of some deviations from 
optimum conditions, concrete, physically realizable solutions were obtained in [14]. These deviations consisted 
of the following: the error coefficients were given beforehand as end values, as a result of which the reproduc- 
tion error, specified by the regular useful signai, was not equal to zero; a graphical approximation to the optimum, 
transcendental, transfer function of the system was introduced by means of the required fraction-rational function 
of the minimum phase type. Another attempt was made in [15], where the optimum solutions for the case of the 
slowly varying (in time) regular signal, were obtained with the help of the minimum, practical, limiting error of 
reproduction criterion. This paper solves the problem of how, starting with the optimum, but usually nonrealizable 
characteristics, obtained in [10], one can determine characteristics which are physically realizable and which 
have the least deviation from the optimum ones. 


1. Statement of the Problem 


Let a useful signal (1) act on the input of the automatic control system in question, during which s(t) ap- 


pears as a stationary, random function having a known autocorrelation function Rss(T) or spectral density G.{w). 


We assume, moreover, that the function of time > af, (t) 1s given, and is chosen on the basis of the most 
k=1 


unfavorable or typical operating conditions of the system, The noise n(t) occurs as a stationary, random func~ 


tion with a given autocorrelation function Rpn(T) or spectral density Gpp (w), and we assume for the sake of sim- 
plicity that the functions s(t) and n(t) are not mutually correlated, 


Designating the pulse transfer function of the linear, automatic control system by w(T), we can write the 
response of the system to the input 6;(t) + 6p(t) in the form 


T 
B, =( w (t) [Bs (t — 1) + Ba (t de, 
0 
where w(r) =Owhen0 =r 2 T. 
Let us designate by w4(r) the ideal pulse transfer function of the automatic control system (ACS). Then, 
in the absence of noise, the response of the system to the input signal 6;(t) will have the form 


Bot wy (2) — 2) de, (6) 


—oo 


The problem ts reduced to the finding of a function w(t) for which minimum dispersion occurs in the 
control errors: 
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83 (t) = E (t) — Bos 
The total control error will be of the form 


T 
8, (t) w(t) [s(¢—*) + n(t—2)] de — \ w, 8 (t—*) de + 
0 —o 
T 


(\ w(t) fe (t de fe (1) 


The last summand in Eq. (7) determines the systematic control error. In [10] the requirement for the sys- 


tematic error to be equal to zero {is introduced, which, on the basis of (7), will be expressed mathematically in 
the form 


T 
0 


—oo 


Following [14] we shall suppose that the systematic control error 5g,(t) isnot equal to zero,and, therefore, 
Eq. (8) is not satisfied. In this case we can write the following set of equations: 


T 
\ (+) fn (t—2) dt w(t) fn (t—2) dt = An(t), 


whence the systematic error is 


n 
= ayA, (t). (10) 
-1 

Because of the technical requirements which are imposed on the control system, the maximum permissible 
values of the functions A(t) (k = 1, . .., m) can be often chosen in such a way that the inequality 158.lmax <€ 
is satisfied, where € is the maximum permissible value of the systematic control error. Subsequently, the pro- 
blem is reduced to finding, for the given functions R,,(T), Rnn(T), Gss(w), Gnn(w) and the observation time T, 
the pulse transfer function w(r) which satisfies the set of equations (9) and which guarantees a minimum in the 
dispersion of the random error: 


T 


0 


T 

\ wy — 2) de — | fx (t—2) de = Ax (0), (9) 
0 
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For simplicity we shall assume that the mathematical expectation of the functions s(t) and n(t) are iden- 
tically equal to zero. 


2. Determination of the Optimum Control Characteristics 


The solution of the variation on the problem in [10] formulated above reduces to the necessity for the 
following integral equation to be satisfied: 


T n 
(+) [Ras (t — 2) + Run (t = + 


where , is a Lagrange multiplier. 
The solution of Eq. (12) has the following form: 


n 2m 
w(t) = >) Aut! + Bee™ + R(p) x 
k=0 


k=1 
x (jn) R(— fo) + 


O<t<T), 


w(t)=0 O>t>7), 


where Ao, Ay,..., An» By, Bg, . . Bym, Cy, Cg... -» Cz Dg, Dg, . . Dz-m are constants to be deter- 
mined, 6(t) is the delta function, 5(Y t) is the vth derivative of the delta function. 


Gas + Gun (0?) = (je) (14) 


2m is the power of the major term of A(w*), 21 is the power of the major term of B(w*), p==d/dt, W; (p) = 


= \ w,;(t)e-dt 1s the ideal control operator and 


s € m 
Yo + +°**YmP 
= = (15) 
he 


The polynomials Q(p) and R(p) do not have any zeros in the right half-plane of the variable p. 


Substituting w(t) from (13) into (12) and treating the resulting expression as an identity, we obtain 21 
linear, homogeneous equations, which contain the unknown constants A;, By, Cy, and Dy, Substituting w(t) from 
(13) into (9) with k = 0,1,..., m, we obtain anothersetofn +1 eqs, The solution of the obtained 2/ +n+1 
equations yields the values of all the constant coefficients. 


The optimum transfer function of the equation of a closed system is, as is well known, of the following 
form: 


| 

(12) | 
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K (p) = \ w (t) e—Pldt. 


Substituting the expression for the pulse transfer function w(t) in (13) into (16) we obtain 


Tn 2m 2m 
K (p) =\ + eon 


T 
4 G,,(@ 
+ R(p) W, (jo) R (— jo) edo + 
0 


4 Cy + Cop Crimp! 4 (Dy + Dap + 


The solution of Eq. (17) shows that the transfer function of a closed control system depends on p and e" PT, 
This reduces to the fact that the frequency characteristics corresponding to the system contain undamped fre- 
quency components, A graphical method of approximating the amplitude and phase frequency characteristics 
of the optimum system with the help of the corresponding characteristics of desired, standard, control systems 
was used in paper [14] to obtain a physically-realizable solution of the problem. 


In the present work an analytical method of approximating directly the transfer function K (p, e~PT) and 
not the frequency characteristics is proposed. To obtain the function K(p) in a fraction-rational form, we approxi- 
mate e P by means of Pade's fraction-rational function [16]. In this case we can write that 


lim 


1 
(p+v)+ 0 Guy () 


var v(v—1)2* 
v(v—1)---2x1x2" 
+ v)(u + v—1)---(v+ 1) 


(4 + v) (u + v—1)-+-(v+ ul 


=1— 


(20) 


Taking x = —pT and p =v so that the moduli of the left-hand and right-hand sides of Eq. (18) be equal 
for any value of yu , we obtain the following approximate relationships; 


forp =v =3 


0 
where 
(19) 
and 
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~ = 


(22) 


The approximate relationships in (21) and (22) are quite adequate for a whole row of practical problems. 
To evaluate these approximations let us examine the phase frequency characteristics of the left and right sides 
of Eqs. (21) and (22). The phase frequency characteristic for the function e~P™ has the form: 
(23) 
for the right side of (21) 


and for the right side of (22) 


() = arctan 


The phase frequency characteristics represented by (23), (24), and (25) are shown in the figure. As is ap- 
parent from the figure, these characteristics coincide very well in the range of frequencies 0 < wT < 6, which 
is the most effective in many control systems. Thus, 
0 2 4 6 8 2-wT in tracking systems, for example, the indicated range 
includes, to a first approximation, all the working fre- 
quencies up to the cutoff frequency. 


The approximate characteristics are different 
from the optimum ones only to the extent that, for fre- 
quencies greater than the cutoff frequency w,, they 
have a smooth shape and do not contain frequency com- 
— ponents brought about by the discontinuous nature of 

netxd the pulse transfer function (13), It {s very important 


= 


to note that the described analytical method of approxi- 
mation, as opposed to the graphical method in [14], 
yields good coincidence between the realizable charac- 
teristics and the optimum ones, within the working range 
of frequencies, irrespective of whether the systems being 
approximated are minimum phase types or not. The final operation in the determination of the optimum charac- 
teristic of a linear control system 1s the substitution of e~PT from Eq. (21) or (22) into Eq. (17), from which we 
obtain K(p) in the form of a fraction-rational function. a. ee 
sating devices may be utilized with a check of the system's stability. 


3. Examples of the Method's Application 


1. To determine the optimum transfer function of a tracking system, given Gpp(w) = ¢ = 0.612°10° 
sec, Ges(w) = 0, Wy(p) = 1, T = 1 sec and given the error coefficients cy = 0.005 sec and cg = 0.0016 sec*, which 
satisfy the requirements of dynamic accuracy. 


for p =v =4 
T (co) = arcean 
0) 
ual 
21) 
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In accordance with Eq. (13) the pulse transfer function w(t) has the following form: 


w(t) = Ay + Ajt + (O0<t<T), (26) 


w(t)=0 (O>t>T7). 


Equations (9) will, in this case, be written as 


| =i, 

0 
| (27) 
0 


T 
(t)dt = — cy. 
0 
Substitution of (26) into (27) yields the following system of equations; 


+ + =a, (28) 


A= — (29) 


Using Eq. (16) and taking into account (26), (29) and‘(21) we obtain the expression for the transfer func- 
tion of a closed tracking system: 


(30) 


The transfer function of the open system {s 


whence 
= 
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For the system to be stable and to have minimum phase angle when in the closed or open condition, it is 


necessary to impose the following requirements on the error coefficients; cy < ; and ~—Tcy< cg< z —Te,. 


The mean value of the square of the error, determined by the interferences having a spectral density 
Gnn(w) = @, will have the form 


2 
Substituting numerical values into (31) we obtain 
0.0967 p? +- 0.495p 4 


In accordance with (32) we have V %., = 0.093° . The systematic tracking error will be determined 
from the equation 


(t) =e (34) 


Then, in accordance with (33), the maximum value of the systematic error will be (5g.(t))max = 0.192”. 
Examination of Eq. (33) shows that the tracking system is conditionally stable and is astatic in the first order, It 
should be noted that if the first three coefficients of the tracking system's error are given, then it is necessary to 
use the approximation of the function ePT in accordance with Eq, (22). In this case the optimum transfer func- 
tion of a closed tracking system has the form: 


Let and \— 


K (p)= 


and in the choice of the error coefficients the following conditions must be met; 


2. To determine the optimum transfer function of a differentiator, given that Ran (T) = eval rl, Gnon(w) = 


- & = 20 sec™*, G,,(w) = 0, Wy(p) = p, T = 10 sec, By(t) = a + and the systematic differentiation 


error must be equal to zero. According to Eq. (13) the pulse transfer function ts 


w(t) = Ag+ Ayt + C18 (0) + O<t<7). (36) 
w(t) =0 (0>t>T7). 


| 
| 


Equations (9), in this case, have the form: 


T T 
=o, (37) 
0 0 


Substituting (36) into (37) and (12), and treating the last equation as an identity yields the following set 
of equations; 


AT 4+C,+Di=0, + =—1, 


(38) 
Aga — A; — aC, = 0, Aga + A; (aT + 1) — a*D, = 0, 
whence 

+ 6aT + 12 T (a*T? 6aT 12) 

(39) 
6 (aT + 2) 6 (aT + 2) 
T (a*T* + 6aT + 12) ' T (a®7* + 6aT + 12) 


Making use of Eqs. (16), (21), (36) and (39) we obtain the optimum transfer function of the differentiator 
in the form 


+12) 
+ 6aT 12) 10 + + 12) 


K(p)=p (40) 


Equation (40) can be rewritten, after a simple transformation, in the form 


12(aT +2) 1202 4. 120 


T (@*T* + + 12) (41) 
(7+ 58) j (r+ 4; 


whence it follows that the control system is stable. After substitution of numerical values of the parameters, 
Eq. (41) takes the form 


P (p+ 0.3675 0.3505) (pF 0.3675 +7 0.3505) 


The response of the differentiator to a unit step function of speed 8; (p) = 1/p* has, according to Eq. (42), 
the following form: 


Bo (t) = 16,64 1-465! 5.64 ¢—0-3675! (29s 0.3505 + 0.0322 sin 0.3505 t). (43) 


If the interference at the input is determined by a function of spectral density as Gnn(w) = e*, then the 
optimum function of the differentiator can be represented in the form 


K = 
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and the mean square differentiation error will be determined by the expression 


SUMMARY 


1. An approximation of the transcendental function e~PT was used in the work to obtain, with the help 
of Pade's fraction-rational function, the optimum transfer functions of automatic control systems, determined in 
[10, 14], directly in a form suitable for engineering application. 


2. The cited examples point out the practical applicability of the proposed method of approximation to 
the solution of problems in the optimization of automatic control systems, in the presence of interference and 
with given requirements for dynamic accuracy. 
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THE INFLUENCE OF LINEAR ZONES AND REGIONS OF SATURATION ON 
THE DYNAMICS OF TWO-STAGE SERVOMECHANISMS 


N. S. Gorskaia 


(Moscow) 


The dynamics of the free motion of a two-stage servo- 
mechanism are investigated, where the first stage has a relay 
controlling element whose characteristic contains a loop and a 
dead region, and the second stage contains a coutrolling element 
with a linear zone and with saturation regions. The investiga- 
tion is carried out by means of the method of point tansforma- 
tions [1-3]. 


A complete study is made of the nonlinear problem in- 
volving the influence on the free oscillations of the servome- 
chanism of the second stage characteristic with the linear zone 
and the regions of saturation. Analytic expressions are given 
for the critical values of the basic servomechanism parameters, 


I. Posing of the Problem ‘ 


Two-stage servomechanisms find wide application in automatic regulation and control systems, and in 
follower systems, for significant amplification of the controlling signal sent from the sensing and transforming 
devices to the regulating organs, Such servomechanisms are used for amplifying power and torques. 


| supply, supply 
zi | | 
1st stage amplifier 2nd stage amplifier 


Fig. 1. 


The fastest acting servomechanism is obtaned when its controlling element acts with relay, or nearly relay, 
characteristics. For a given structure, stability may be guaranteed either by introducing dead zones, or by intro- 
ducing linear portions into the “yes-no” characteristic of the controlling element, 


The introduction of dead zones leads to a decrease in the accuracy of the servomechanism; the introduc- 
tion of linear portions in the relay characteristic allows accuracy to be maintained. Since the size of the linear 
regions, when the aim is to increase the quality of servomechanism performance, is chosen to be small, when 
the computations are made for such a servomechanism the compound characteristic of the controlling element 
is ordinarily replaced by a relay characteristic. 


In the majority of published works, servomechanisms 
controlled by only one relay are considered, and only 
certain works [4, 5] are dedicated to two-stage servo- 
mechanisms, In the works cited, the second stage approxi- 
mation of a servomotor is considered as; in [4], a relay 
characteristic with a dead zone, in [5], a purely linear 
characteristic. 


Here, more typical servomotor characteristics will 
be considered, characteristics described by uneven non- 
linear functions with linear zones and saturation zones, 


Fig. 2. Schematic of a two-stage pneumo-hydraulic 
servomechanism: 1) the controlling element (the noz- 
sle-valve type), 2) membrane executive mechanism, The problem of the present work, from the technological 


hyd point of view, consists of obtaining recommendataions 
it + for the rational construction of two-stage servorechanisms, 


the block schematic of which is given in Fig. 1. From 

the mathematical point of view, the problem comes 
down to the investigation of systems with relay characteristics containing loops, and with nonlinear characteristics 
with linear and saturation zones. This allows a comparison to be made of the influence of the linear zone p; 
and the dead zone pq on the system's stability and auto-oscillation. The investigation is carried out by the 
method of point transformations. 


As an example of the type of servomechanism considered, one may adduce the pneumo-hydraulic servo- 
mechanism (Fig. 2) [4-6]. 
Il. Equations of Motion 


The equations describing the free oscillations of the two-stage servomechanism, with the assumptions of 
[4, 5], but taking into account the saturation zone and the linear portions in the nonlinear characteristic of the 
controlling element of the second amplification stage, may be written in the following way: 


Tye +p = (2), = ;(p), (1a) 


where * p is the relative magnitude of the output coordinate of the first amplification stage (relative pressure 
variation in the housing of the membrane mechanism), x is the relative magnitude of the output coordinate of 
the servomechanism (relative displacement of the servomotor shaft), T; is the time for establishing a pulse in 
the first amplifying stage and T, is the response time of the servomotor. 


The nonlinear functions, (x) and 3(p), are represented as: 


1 for if O(a) =4, 
(2) ={ 0 for =0, (1b) 
for tf = —1 


(where x, is the initial condition at time t = +0), 


1 for P> Pp 
p for (1c) 
—1 for P< — Py, 


* Since we are considering the free oscillations of the servomechanism, the input signal, the rotation of the mani- 
fold with nozzle with respect to the valve, equals zero, and the opening of the pneumo-relay valve depends only 
on the action of the feedback path (5). 
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ait 
pressure 


where A, 4 and p, are relative quantities for the slack in, respectively, the feedback path, the dead zone and 
the linear region of the relay. 


III. Salient Features of (p ,x) Phase Space 


The (p , x) phase space (plane) in the problem under consideration differs from the (p , x) phase plane con- 
sidered in [4, 5]. 


AL 


N 


Letting 
wa 


Fig. 3. Three-sheeted (p , x) phase plane and limiting cycle for the 


In the given problem, the (p, x) space, in correspondence with the characteristics of ©,(x), is a plane com- 
posed of three sheets lying one on top of another, each of which, moreover, in correspondence with the charac- 
teristics of 3( p), is divided into three regions by the auxiliary lines, p = +p, (Fig. 3). 


The equations of the integral curves in the various regions are written as follows; 


in region 1 
tin — ug (2) 
where Ts 
in region 5 
(1 + py) — In (4 + = + (3) 
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oxi- 
a 
| 3 
case i—p,e — + 
ni- 
y 


t{[—p, — In (4 — — — In (1 + = — + 


(1 + peg )— In (1 — = — 


Here, the numbers in the subscripts refer to the 
region number, and the letters "i" and "f* refer to the 
initial and final values of the coordinates in the region. 


If there is a periodic oscillation in the system, 
Paf = ~Pxi- Then, based on Eqs. (2)-(5), we have 


Pri = * = — Pt = 


\2 
e (1+, 
1+Pi¢ 


We obtain from (6) the equation expressing the 
condition for existence of periodic oscillations in the 


system: 


A 
md 


1—pfe* = 


If Eq. (7) has real and positive roots, p sf, then 
there exist periodic oscillations in the system. We 


Fig. 4. Diagrams of point transformations. 


* = 


A+ Uo (pif). 


The curves of the functions u(p sf) and ug(p sf), represented in Fig. 4a, have only one point of intersection 
for 


"le 


| in region 3 
(4) 
in region 4 
uy 
4, (6) 
b 
. a 
~ 
ty i+ f (7) 
\ 
\ 
|_| 
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If the relationship of the system parameters is such that 
1 =e (1+ (10) 


then p47 = Pg, and a limiting cycle {s set up on the phase plane, represented in Fig. 3, and auto-oscillation 
arises in the system. 


1. Weconsider the (p , x) phase space for p; < ae (Fig. 5), containing the following integral curves; 


in region 1 


4. 
’ 


in region 3 
= pp — In (1 — — Im (1 + = — teat 
in region 4 
= {In (1 + pet) — In (1 — = — — 


Periodic oscillations in the system. In the cases where periodic oscillations are present 


A 
Pag = — Pri and = — 


Based on Eqs. (11)-(14) we have 


(1 + prog) = 


We introduce the notation 


4 
+ pag) = = Us (Pr) (18) 


e 


(uy (P gf) is a linear function, ug(p gf) is the sum of two exponential functions, ugy(p gf) and uge(p 2f)}, and we shall 
seek a graphic solution of Eq. (17). 


The curves uy(p gf) and ug(p g¢) intersect the axis of ordinates (cf. Fig. 4b) in the points 


{tts (Pag {ta = [1 ] > 0- 


With this, { ug( Pat) > for > 0, A>0,p,> 0. 


From an analysis of the functions it follows that the curves u;(p gf) and ug(p gf) intersect in only one point. 
For each ratio of the servomechanism parameters higher than the critical ratio, uy (p9f) and ug (P9¢) intersect only once 


in region 2 
— = + iE: (12) 
(13) 
(14) 
py e°2f 
(16a) 
e 
or 
j 
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(P af = Pap), 1. €., for the case p, < “ » periodic oscillations can arise in the system only with a definite ampli- 
tude Py. m4 


ratio of the It follows from system Eq. (16a) that pe = 0 forpy= pij = * 
where p > p; since a<—. 


This indicates that the boundary of the dead zone in the first relay is reached at the moment of switching 
in of the second controlling element. 


Thus, a necessary condition for the existence of periodic oscillations is the holding of the inequality 
Pry > (Or > 
Substituting p gf = 0 in Eq. (17), we obtain the following critical ratio of the parameters; 


ake 


T=i—pe* 


ALLL 


2 


Fig. 6. Diagrams of the point transformations of the Fig. 7. Diagrams of the point transformations of the 
half lines N and N* to the segment M"M" of line M for _half-lines N and N° to the segment M"M™ of line M: 
the case < 2uq/r: a) for tg < b) for =a) for the case p, < (for #4 > >) for 
Hd 


the case 1 me ¢ 14 + 


ty) 


The value pt g = #4 or yy» for which Equality (20) holds, is the critical value. With this 


Uy Uy, Ug ly 
if 
2 
| Bil 
4 / 
4 
4 
hall 
nt. 
ly once 
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The ratio of the parameters determining the boundary of the qualitatively different cases of servomechanism 
motion has the form: 


A 


hq 
e <i—pe* (21) 


A. The construction of the diagrams of the point transformation for tg < Hdcry (for the case p, < 


2 
< ts considered in Appendix I 


The curves of the functions ug(p 31), Ug(P sf). Us(P 21) and ug(p yi) (Fig. 6a) have one point of intersec-ion 
Py, and the resulting limiting cycle will be stable. 


_ The amplitude of auto-oscillation, depending on the output coordinate (x) of the servomechanism is deter- 
mined from Eq. (13), by substituting in it xs = X , Per = 0, and from system Eq. (16a), by substituting in it py = 


= Dy: 
—Ing 


——— Sheet Z 


Sheet NW sheet I 


4 
| 
| 
| 
| 


ro 


B, Construction of the diagrams of the point transformation for # 4 =  g cy yy» 1. €., when the ratio of the 
parameters is determined from the equation 


2hq 


4 


is considered in Appendix II. 


It follows from the diagrams of the point transformations (Fig. 6b) that curves us(p 4) and ug(p y) intersect 
for px = pn. 1. e., for the limiting value of py for which periodic oscillations of the system are possible (for 
2 


In this case, the limiting cycle formed in phase space (Fig. 5), as follows from the diagram of the point 
transformation (Fig. 6b), is semistable. 


This latter means that the least variation of the parameter ratio to one side or another (of the critical ratio) 
will lead either to the destruction of the periodic motion and the damping of oscillation or, conversely, to the 
formation of a stable limiting cycle. 


Sheet 7 


Fig. 9. 


C. The construction of the diagrams of the point transformations for > gory when the ratio 
of the parameters is determined by the equation 


"t>1—pe* 


P 
Sheet 
| 
is considered in Appendix III. 
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Figure 7a shows the relative placement of the curves us(p 34), Ug(P af), Us(P si) and ug(p yf), allowing the 
conclusion to be drawn that, for 4 > 4 .. 17» periodic oscillations do not arise in the system. 


For #a> Maern the system, after an arbitrary initial deviation, moves to an equilibrium position. 


D. The construction of the diagrams of the point transformations in the case when the ratio of the system 
parameters is determined by Eq. (10) is investigated in Appendix IV. 


Figure 7b shows the relative placement of the curves us(p 4j), U4(P 11), Us(P x1) and ug(p y), allowing the 
conclusion to be drawn that, if the system parameters are connected by Relationship (10), then a stable limiting 
cycle is formed on the phase plane, and auto-oscillation is established in the system. 

2 

2. > It would seem (from a comparison with the case 

Qu 
Pi< a. ) that, with this ratio of the ee Pe parameters, the least value of the initial deviation of 


the ais ME p for which auto-oscillation is possible, must satisfy the inequality p it <p, . However, as investi- 
gation showed [5] 


2ud 
(23) 


2 
The semi-stable limiting cycle for the case pz >as is shown in Fig. 8. 


it follows from [5] that only one limiting cycle can exist in the given case. 


2u 
For p; > —, the parameter ratio determining the boundary of qualitatively different cases of servo- 


mechanism motion has the form: 
(24) 


3. We now consider the case when the ratio of the servomechanism parameters is determined by the equality 
Py = ate, The limiting cycle for these parameter values, as represented on Fig. 9, is a “bounding” one with 
respect to the limiting cycles for Py < <td and p; > el It is easy to convince oneself of this by consider- 
ing the poisition of the limiting cycles mentioned on the (p , x) phase plane. 


The parameter ratio determining the boundary of qualitatively different cases of servomechanism motionfor 


IV. System Auto~-Oscillation 


The investigation of the dynamics of a two-stage relay-type servomechanism containing relay functions 
with loops and uneven nonlinear functions with linear and saturation zones showed that it is possible to have stable 


2ud A + 2der 
© > for < Hderl defined by ihe equation exp 
Hderl 
or if py < —— for < ts defined by the equation exp (- 


T 


If in the servomechanism the magnitude of the linear zone p; > then for < Mg cry, & limiting 


cycle is formed on the phase plane inside the linear zone. The amplitude of auto-oscillation is Py < p;. If 
however then the amplitude of auto-oscillation is Pyj > pl. 
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2 

pp =— has the same form as for py > ——*: 

(25) 
=1— p; exp 1 


2 
If the width of the servomechanism’s linear zone {s P,< =e then, for #q < Hg crys @ limiting cycle 


is not formed inside the linear zone on the phase plane, since the previous trajectory, with initial deviation 
Pili < pz, will tend to a quiescent segment, and auto-oscillation may arise only outisde the linear zone, If 
now Hg > Udery there will be no limiting cycle even external to the linear zone (for p yj > P;) and all tra- 
jectories will be deflected to the quiescent segment for any initial deviation, p yf. 


And, finally, if tg < Mg oryye @ limiting cycle will be formed leaving the linear zone. 


SUMMARY 


1. A two-stage servomechanism with a relay characteristic with loops in the first amplifying stage an un- 


even nonlinear characteristic with linear and saturation zones in the second amplifying stage has two critical 
parameter ratios. 


I- A Pa In the case when the width of the linear zone ts 
2 


Py ed , the critical parameter ratio has the same 


form as in [5}. If the width of the linear zone is p; < 


2 
< = , then the critical ratio has the form of (20). 


2, The presence in the second stage of a nonlinear 
characteristic with saturation leads to the dynamics of 
such a servomechanism being essentially different from 
the dynamics of a servomechanism with a purely linear 
controlling element in the second stage [5],or one with 
a relay with dead zones [4], 


3. If, for a given steepness of the linear portion 
of the controlling element's characteristic, the system 
possesses auto-oscillation leaving this zone, then the 
amplitude of the auto-oscillation will be less than for 
a purely linear characteristic. This is explained by the 
decrease in the steepness of the characteristic with in- 
crease of the coordinate p. 


ality 


ma The introduction of saturation regions in the charac- 

teristic of the controlling element of the second stage 

leads to a decrease in the critical value of the dead zone, 

if the linear zone is selected by the condition that pj < 
2ud 


003 006 209 


| 

| 

| 

T 

auto- oscillation 

ifor 

1 


4. The critical ratio of parameters in the case 
Fig. 10. 2ud 
Pee EF , defined by Ratio (20), has the form: 


There was obtained in [4] the critical parameter ratio (Eq. (3.13)) for a similar structure of a servomechanism, 
T this one, however, having a relay controlling element in the second stage with a dead zone py. This ratio is 
4 


defined by the equality p , = 


e*—1 
ting It follows from the dependence of the critical values of p, and p, on Wy= for A* = © = 0,001, 


423 


as given in Fig. 10, that the linear zone narrows the region of auto-oscillation for values of the first-stage dead 
zone of tg > 0.04. 


For us, < 0,04, the second-stage linear zone broadens the region of auto-oscillation, in comparison with 
the auto-oscillation region for the dead zone p 4. 


The present work was carried out under the direction of V. V. Petrov. 


Appendix I 


2 
We consider the diagrams of the point transformation for 4 < H dcr 1 (for the case pz < st). We 
show that here the equilibrium position is stable when the initial deviation does not exceed a definite value 
(for example, py > forx = — 


We shall also show that when periodic motions occur they will be stable, i.e., that for wy < decry the 
initial deviation exceeds p yj. For this, we consider the mapping of a half-line onto a half-line. 


For the given case, it is expedient to study the transformation of the half-line N onto the line M, and the 
half-line N’onto the same half-line (Fig. 5). 


The half-lines N and N° are parallel to the axis of ordinates and start from the points with abscissas, re- 
spectively, -(a- $ and wg+ . The line M 1s parallel to the axis of abscissas and starts from the point 
with ordinate p;. On the basis of Eqs. (3)-(6) we have 


[ ] 


i+po 


Imposing Condition (15) on Eqs. (26) we obtain 


= (Pj) = Us (P,; 


A 1— py 


‘If Condition (21) is met, the functions ugs(p sj) and u,(p yj) are the expressions for the functions correspond- 
ing to the point transformations of half-line N and half-line N* to the segment M'M" of line M. 


We now consider the mutual placement of the curves of ug(p 3) and u,(p yi) (Fig. 6a). 


2 
It was previously proved that, for p; < —=< , perlodie oscillations could exist in the system only in case 


*t and the parameter ratio satisfied Condition (21). 


For py = 1, we get 


| 

(26) 

A i 
4m 


A 4 
Us (Py) = t+ U4 = 00, 


2 
from which we conclude that if Inequality (21) holds and for p, < a. 


A 
> < {us og < {Me (ip 


When the {nitial deviation of coordinate p equals ply =e,¢* , the function us(p qf) attains its limiting 


A A 
The function ug(p attains its limiting value, equal tou 4 + = 


for 
Py * (+e). 


It follows from the orientation of the curves of us(py)and uy (py) that 


) Us 
In fact, if Inequality (21) holds, 


el, = >i-e * T>,e* = Py 


*d 
eye? * 


and if the transcendental function e?? is expanded in a series, limited just to the first two terms, we obtain only 
a stronger inequality, i. e., we may write 


m a 


—29 
Thus, an analysis of the functions us(p yj) and u,(p 41) shows that the curves us(p qj) amd ug (p do not inter- 
sect in the interval pf < py < 1 if the parameter ratio satisfies Condition (21). 


The point transformations of the half-lines N and N' to the segment M"M™ of line M will be carried out 
by means of functions us (p 44) and ug(p yj) which, based on (11)-(15), may be written in the form 


= = = 4s 


A 
+ — 1) = Paf 
= Ue 


or 
A 

(28) | 
| 
t+ ° 
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We now determine the placement of curves ug(p 4;) and ug(p 4i), as described by Eqs. (28). 


Since it was previously proved that only one limiting cycle can exist in the system, curves us (p 4j) and 
Ug (P 4) can intersect in only one point. 


Moreover, it was proved that periodic oscillation of the system (for P,< - can exist only in case 
2p, 


Pyj > the parameter ratio satisfies Inequality (21). 
Functions ug(p 4) and ug (p yj) are identically equal for py = py. 


Functions u,(p and ug(p yj) are identically equal for py 


With an initial deviation of = Py] = * the function us (p 1) equals + — 


At = Pz the function ug(p yj) attains its limiting value, (us —(Ha— +) 


The function ug (p 34) attains the same value, for 


A 


It follows from Eq. (29) that 

since, if we take the expression (1 + p j ~ —~~) as an approximation to the transcendental function ¢ . 

by means of a series carried to only two terms, then 


Then if Inequality (21) holds, 


Ill lll 
Pig > andes; > 


Thus, the orientation of the curves us (p yj) and ug (p yj) allows us to conclude that the curves have a point 
of intersection for p sj = Pi, and that the resulting limiting cycle will be stable. 


Appendix II 


We now consider the diagram of the point transformations for yt 4 = # 4 or 11 (Fig. 6b), 1. €., when the 
parameter ratio is defined by Eq. (21"). 


The mapping of the half-lines N and N’ into the segment M'M” of line M is carried out by means of the 
functions ug(p 4) and u,(p 4), described by the system of Eqs. (27). 


For p 3{ = 1, in the case when the Equality (21°) holds, and for p; < —*, we have 


A 
gt < {us < {ua (Pri Yo, i-1 = 


With an initial deviation of py, equal to pi; =p, ¢* , the function us(p 4) attains its limiting value, 
ward. The function u,(p 3) attains its limiting value, equal to (u q + S ), for P ai, equal to Pa s 
= 
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a 


From the placement of curves us(p 4) and u,(p 4) it follows that 


> 


since, if Equality (21") holds, 


and if the function e”! is expanded in a series limited to only two terms, we obtain only a stronger inequality,t.e., 


4 
— 
The mapping of half-lines N and N* into segment M"M™ of line M is carried out by means of functions 
Us (P 4i) and ug(p 1), described by the system of Eqs. (28). 


The functions ug(p yj) and (p xf) are identically equal for = pi. The functions ug (p and ug (p 
are identically equal for py = py. 2g 


For an initial deviation of we have 


4 


From (28) we determine the value of ug(p si) for pi = Py: 


4 
i—p,e* 


Whence, if we expand e~P 2f in a series limited just to two terms, we get 


{ue <" {us (ip, 


Appendix III 
We consider the he > 1. when the parameter ratio 
Hd 
is defined by Eq. (21")(Fig. 7a),. Then, for py = 1 and p,; < + we have 


4 
+ < {us (Oy Notes < {ug (Py Nou da = 


For and the same conditions, we find 


_4 
or 
4 
>t—e 
Pot 
nt l 
4 
For Equality (21") to hold, p af = 0, whence {us  and,consequently, 
iu 
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A 
{ug 


The function u,(p 44) attains its limiting value, equal to pq + ; , for 


- 
* 


If Inequality (21") holds, 


Yd 


- 


If the transcendental function e?! is expanded in a series in which only the first two terms are retained, 
the sign of the inequality is reversed: 


ee 


We now consider the placement of curves us(p yi) and ug(p yj), described by Eqs. (28). The functions 
gi) and us(p yy) are identically equal for pj = The functions u,(p and ug(p are identically equal 
for = 
For =F e* wehave 


| 


A 
{us = +3 Te), 


Whence we conclude that, if Condition (21") is satisfied, 


{us (64; > {us 


Appendix IV 


In analogous fashion we consider the diagrams of the point transformations when the parameter ratio is 
defined by Eq. (10) (Fig. 7b). 


We analyze the placement of the curves us(p 4) and u,(p 44), described by the system of Eqs. (27). 
For py = 1, if Equality (10) holds, 


A 
< {4s Oy < {Us = 


For Py we have 


A 
i—e * oft 
or 
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4 
{ug (ey Nog == {uy (244)! = + 


Pu 


We consider the functions us (p yj) and ug (p yj) which are equal to each other, and to the functions us(p 4) 
and (p x1), for px = ply. 


For pj = ;, the function ug(p sj) attains its limiting value: 


A\ 


The function ug (p sj) attains the same value, — (a-$ , for the value Pt = pill , which 1s determined 
from Eq. (29), from whence it follows that 


_4 
If the function e”? is expanded in a series and only the first two terms are retained, we obtain 


4 
elt <i—e * 


or, if Equality (10) holds, 


| 
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ANALYSIS OF FREE OSCILLATIONS ABOUT ITS CENTER OF GRAVITY OF 


A NEUTRAL PLANE WITHOUT DAMPING OF ITS OWN AND WITH A 
RELAY AUTOPILOT*®* 


V. Iu. Rutkovskii 


(Moscow) 


A relay system, the linear part of which is described by 
the simplest degenerate third-order equation, is considered, The 
control system's dynamics are analyzed by means of the method 
of point transformation of surfaces, which enables one to show all 
the possible modes of movement in it, including the sliding modes. 


The analysis yields the equations of surfaces determining 
the region of attraction (by the initial conditions), equilibrium 
state, and the stable limiting cycle. These equations may serve 
as criteria in the choice of the regulator's parameters, under given 
restrictions, from a coordinate and two of its derivatives. 


1.-Statement of the Problem 


The present work is devoted to the study, by means of the method of point transformation, of the simplest 
relay system, the linear part of which is described by a degenerate third-order equation. A similar but more 
general problem was studied by A. A. Andronov and N. N, Bautin [2, 3], B. V. Bulgakov [4], A. M. Letov [5], 

V. A. Kotel'nikov [6], P. V. Bromberg [7], I. Flugge-Lotz [8], etc. However, results of the study made on the 
problem being considered cannot be obtained in a particular case from the above papers. In fact in [2] and [3], 
the authors make use of such generalized coordinates, intrinsic parameters, and dimensionless time, which in 
our case (M = 0) all become zero. 


It should be noted that in the case where the solution is carried out in usual coordinates, it is possible, in 
principle, to obtain results in the study of a system, the linear part of which is described by a degenerate equa- 
tion, as a special case of the study of a system described by a complete equation. For this it is necessary to carry 
out a transition in the limit. However, this problem usually turns out to be more complicated in comparison to 
an independent study of the system described by the degenerate equation. 


The results of [8] cannot be applied to the solution of the given problem not only in view of the above- 
mentioned complications in applying them, but also because of their incompleteness. 


It follows from the papers of the other authors that the system under study is mostly unstable. However, 
the limits of the system's stability, determined by the initial conditions, have not been found, and there is, there- 
fore, no criterion for the choice of the regulator's parameters, It should be pointed out as well that it does not 
appear to be possible to find the regions of stability from the initial conditions using the method of the small 
parameter and the frequency method, while the solution of the given problem by the method of harmonic balance 
yields principally incorrect results. Thus in the study of the system by the indicated method, we are led to the 
conclusion that the blind zone in a servomotor leads to oscillations. However, as will be shown further on, the 
blind zone merely decreases the region of stability determined by the initial conditions. 
* The work was carried out under the direction of B. N. Petrov. 


Thus it would be expedient to solve the given problem by the method of point transformation which would 
permit, in general, the derivation of the equations of surfaces determining the region of attraction (by the initial 
conditions), equilibrium state, and the stable limiting cycle. The latter permits the choice of the regulator para- 
meters in such a way that the system {s, under given limiting conditions, efficient. 


Solution of the problem {s based on [2, 3], These were the first papersto appear in which the space problem 
was rigorously solved and regions of stability were established for it in terms of the system's parameters, The 
theory of the point transformation of surfaces is set forth in them as well. 


2. Equations of Motion 


Let us examine the free oscillations about its center of gravity, from its angle of course, of a neutral plane 
having negligible damping of its own, steered by an automatic pilot with a constant speed servomotor. 


The equations of motion of the system under study 


~@—-— + may be written, after the usual simplifying assumptions, 
as 


Tsp = F(x), 
n=kee+To—m, 


(1) 


where ¢ ts the relative change in the plane's angle of 
course, is the relative change in the angle of rotation 
of the steering wheel, n 1s the relative change in the 
valve's coordinate; k, kg T, y are the constant coefficients characterizing the parameters of the system; F(m) 
is the relay function which, in general, is not single-valued. A schematic drawing of the SAR {s shown in Fig.1. 


Let us introduce new variables = 7'gu, and denote 


N= ~. Then Eq. (1) will take on the form: 


= f (0), o=—p—pp—+8 


z= Nb, = — + f (0), 


where z= 
The relay function f(o) is described by the equation 


1 for o> o,— and/ (a) = 4, 
0 for |o|<o, + -5-and (a9) = 0, 
—1 for (0) = —4, 


where og is the region of insensitivity, A is the loop width, og are the initial conditions at each region of re- 
alignment with t = +0, (At every change of equations the time origin is reset to zero). 
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3. The Stabilization of the Plane’s Course by Means of an Automatic Pilot with a 
Relay Servomotor Not Possessing a Blind Zone 


We shall consider z, 5 , and o to be the rectangular coordinates of phase space, and examine the solutions 
of Eqs. (2) under the initial conditions t = 0, z = z»,5 = 5) and og =». At the same time we assume that z(t), 
5(t) and o(t) are continuous at the points of discontinuity of the relay function f(o). 


Fig. 2. Phase space with o, = A=0. 
Plane A (o = 0) divides the phase space (z, 5 , a) into two parts (Fig. 2), In the top f(a) = 1 and in the 
bottom f(a) = —1. The phase coordinates change in time according to the expression: 


B= f(0)t-+%, z= Nf + + 2, 
— Nf (2) — ot (2)) N —(zo-+ + f(a) +29. 


Projections of phase trajectories on the planes o = oy = const are described by the equation 


N28, N 5 


while the projections of the phase trajectories on the vertical planes z = Zp» = const are the cubic parabolas 
BN Zo N 952 4 
— Her? + Fer 


i Y 
+5 


Thus each phase trajectory appears as a line of intersection of two surfaces, of which one {is formed by the 
movement of the vertical axis according to (5), while the other is determined by Eq. (6) if one substitutes into 
it Z» = Z9(z,5 ) in accordance with (5). 


Equations (2) determine the point transformation of the plane o = 0 into itself. To the extent that phase 
space is symmetrical with respect to the origin of the coordinates, it suffices to examine the transformation of 
the left (5 < 0) half-plane o = 0 into the right (6 > 0)half-plane o =0, We shall call this transformation S*, To find S*, we 
proceed in the following manner; let the starting point z»,5, (t = 0) be in the plane o = 0, Let us use T, to 
designate the smallest value of t at which the phase trajectory, having originated at the point (z», 59), once again 
appears on the plane o = 0. Then, using the derived value of 5», we find z, from the last equation in (4): 


20, T=— FT 


After time t = T, the phase trajectory will cross the plane g = 0 at the point z,;, 5g 


4 =%4(6, 7) 


Equations (7) and (8) yield, in parametric form, the transformation £, 


We shall begin the study of the point transformation S* with the construction of parabola Ly (Fig. 2), which 
possesses the property that the phase trajectory crossing it ato = 0,5 =459< 0, will intersect it again at o = 0, 
after time t = Ty. We shall obtain the equation of this parabola by substituting into the second expression of (7) 
or (8) Ty = 2€: 


(9) 


For f(o) =—1 we shall analogously obtain para- 
bola L,: 


(10) 


/ Since for the existence of an ordinary limiting 
7 cycle in the system it is necessary to have 5; = —5o, 
re J such a cycle must pass‘through the points of intersec~ 
Ww tio of parabolas Ly and L, with the 6 -axis 


(11) 


Fig. 3. Determination of the stability of the limiting Ts 
cycle, + V 


A study of the point transformation shows that in the system, a limiting cycle that passes through points 
(11) is possible. Its projection in the plane o = 0 consists of parabolas L, and Ly: 


FE. (12) 


| 

== 
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To determine the type of cycle, let us first note that the phase trajectories lying in the arbitrary cylindrical 
surface, the projection of which on the plane o = 0 follows, for example, parabola L, (Fig. 3), converge to the 
trajectory crossing parabola 1;. In fact we shall find from the third equation in (4) the expression for T;: 


= = 0) 


N 


T,=3 


If Zp = < 0, 5g =5, < 0 (Fig. 3), then T, = 2]5,]; 59 = 5g" < 0, then = 2|5, 1, 
< 2g", 59 = 5, = then > and hence 5g > . can be that 


Let us also show that the obtained limiting cycle 
is an unstable one, For this it suffices to examine the 


Projection of 
unstable limiting motion of the representative (reference) points which, 


XO: ly cycle in the first instant of time, are located on the switching 


mode 


plane o = 0. Let us take, for example, point g which 
lies outside the area enclosed by the projection of the 
limiting cycle (Fig. 4). After some time t; > 0 the 
reference point travelling along the phase trajectory, 
situated below the plane o = 0, will arrive once again 
on this plane at point d. In this case, in accordance 
with the property just examined, the inequality 5, > 
; > [54] will hold. Thence, the reference point will 
continue its motion along the trajectory, passing above 
the plane o = 0 and, after t, > 0, will once again re- 
een ane turn to plane g = 0 at point e, where 5; > | 5,[, etc. 
: In this way the reference point being considered recedes 
from the origin of the coordinates. A similar argument 
shows that a reference point which begins its motion on the plane o = 0, within the area enclosed by the projec- 
tion of the limiting cycle, approaches the origin of the coordinates, which proves our assertion. It should also 
be noted that the system enters a state of equilibrium z = 5 = o = 0 in a sliding mode, 


The sliding mode in the system is reflected in the motion of the reference point along the lamina G, of 
plane go = 0 enclosed by the straight lines L, and L,: 


t= (14) 


Equation (14) is obtained from the fact that T,; on lamina G has a negative value for f(o) = 1 as well as 
for {(¢) = —1. To determine the nature of the motion of the reference point along lamina G, we must assume 
that the relay switches over with a certain time delay r. Pre-establishing the system in this manner we arrive 
at the obvious result that, assuming symmetry in the projections of the phase trajectories with respect to the z- 
axis, the reference point moves on the lamina G along the straight lines 


8 = 8) = const. (15) 


Upon reaching either L, or Ly, the reference point moves towards the z-axis along these straight lines if 
55 < 5g < 5g», or moves along the corresponding phase trajectories if 5< 55 or 5g > 5s". Here 5s = —B and 
5s* = +8 are the abscissas of the points of tangency of the straight lines L, and Ly with the parabolas (a) and (a"), 
which are determined by Eqs. (5) (points S and S* in Fig. 4). 


In [9] it is shown, that in relay systems, the linear part of which is described by an equation of not less than 
third order, there can appear not only the simple limiting cycles but also ones of any degree of complexity. How- 
ever, in the problem under consideration, in view of the fact that the equation of the linear part appears as the 
simplest form of degenerate equation, complex oscillations are absent. This follows from the fact that the refer- 
ence point, as it was just shown, taken at some instant of time on the plane o = 0, outside the contour bounded 
by parabolas L, and Ly, with subsequent motion will continuously recede from the origin of the coordinates (6 
incessantly increases), On the other hand, if the reference point on the plane o = 0 is within this contour, then 
subsequent motion will only result in its approaching the origin of the coordinates. For this same reason in the 
system under consideration there cannot be limiting cycles consisting of sliding and nonsliding motions, the pre- 
sence of which in high order relay systems has been proved in [10] and [11]. 


In fact, if the reference point, on leaving the sliding mode lamina Gy, arrives, after one half-cycle of os- 
cillation, on the plane o = 0, outside the region enclosed by the contour formed by L, and Ly, then the system 


further on will break into runaway oscillations, In the opposite case it will settle into a state of equilibrium (in 
a slipping mode), 


Thus the system being studied is mostly unstable. The region of stability is, by the initial condition for 
Gz, enclosed in the space between two cylindrical surfaces M and N, shown in Fig. 2, and bounded above and 
below by the surfaces whose equations, in parametric form, are 


Moreover, the system will be stable if the initial conditions correspond to the space Gs from which the 
reference point will arrive on lamina G, under the conditions || < max (6 ,J5*|),where6&* is the pote of inter- 
section of Ly and L¢. 


The region of stability of the system widens with increase of the stable feedback factor 1/a, the reaction 
coefficient along the derivative o(B ), and with decrease of the servomotor's reset time. For objects with a large 
moment of inertia the region of stability is greater than for low inertia objects, To determine the region of sta~ 
bility along the angle » we must make use of the formula 


Yo = — — — (17) 


Thus if the initial conditions belong to spaces Gz and Gy then the system will arrive at the stable condition 
z= 5 =o=0; if they lie outside these spaces, the system will break into runaway oscillations. 


4. Stabilization of a Plane’s Course by Means of an Automatic Pilot with a Relay 


Servomotor and with a Blind Zone 


Phase space (z, 5, consists of three parts: I) (f(a) = 1), II) (f(o) = 0) and IM) (f(o) = 1). A join- 
ing of trajectories, corresponding to different parts of phase space, occurs on the planes A(o = o¢) and B( o = 
=—o,). In parts I and III the phase coordinates change in time according to (4), the projections of the phase 
trajectories on the planes go = gg = const are determined by Eq. (5), andon the planes z = Z» = const by Eq. (6). 
In part II the phase trajectories lie in the vertical planes 5 = 5, = const and are determined by the equation: 


am + BN’, —24)| + (18) 


In view of the symmetry of phase space, shown in Fig. 5, we shall examine the point transformation of 
plane A not into itself, but the transformation of its left half (6 < 0) into theright half (6> 0) of plane B. 


| 
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| 
| 

435 


Fig. 5. Phase space with o, # 0, A= 0. 


Fig. 6. Diagram of point transformation for the case 
€0,4=0, 


Let us call the transfer of the reference point (z», 
59) from the left half-plane of A to its right half-plane, 
transformation S*. Expressions (7) and (8) give the trans- 
formation S* in parametric form. 


The transfer of point (z;, 5;) from plane A to plane 
B with & > 0 we shall call transformation E*. We shall 
define the transformation E* by means of the parameter 
T, > 0, which is the time required for the point to make 
the indicated transition. We call the coordinates of the 
image point on plane B( 5 > 0) Zs, 5, and write the 
transformation E* in the form 


6, = 8, (0, =4,, 


(19) 
= 22 (5, Ty, = + x. 
For a simple limiting cycle to exist in the system 
it is necessary to have 5, = —5,, and z; = —Z». From 
the first condition we find that 


T, = 2, (20) 


while from the second we have 


Let us now find T, from the condition that the reference point (z3, 54), having begun its motion in plane A 
will, after time T3, arrive on plane B (zs, 53). For this we can make use of the third equation in (4), after setting 


f(o) = 0 in it. Asa result, we shall obtain 
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T, Nt Y, (0). 
Comparing T, obtained from Eqs, (21) and (22), we can find the 4 -coordinate of the amplitude of the limit- 
ing cycle; 
\ 
=V (23) 


The relative position of curves Y;(¢) and Y2(¢) is shown in Fig. 6. Examination of these curves, taking» 
into account the results of the analysis for the case o¢ = 0, shows that there is only one unstable limiting cycle 
in the system. The functions Y,(¢) and Y,(¢) may a 

have only two common points; ¢ = 0 and ¢ = 5,. The 

limiting cycle (Fig. 5) goes from the point (z9, 59 = 54) q 

along the cylindrical surface M, then a switch over f(a) i 
from the value 1 to the value 0 occurs on plane A and 


the image point from then on follows the vertical plane 
| C (6 = —&» = const) up to plane B, The second half 


of the cycle is symmetrical with the one just described 


with respect to the origin of the coordinates. 
A sliding mode, which travels either along the 
switching plane A with 6 < 0 and with -BN6 > z > 
> — BN6 — Va, or along plane B with 6> 0 and 
— < 2<—BN6 + 1/a, is possible in the system. 
=. Figure 7 shows the sliding mode in the system for the 
wane- case when it travels along plane A, As can be seen from 
the diagram, the distribution of the regions of the sliding 
mode when o¢ # 0, is significantly different from the 
plane distribution for the case o ¢ = 0. Moreover, the pre- 
all by determination of the sliding mode at the expense of the 
eter delay in the relay shows that the reference point moves 
nake Fig. 7. Determination of the region of the sliding along inclined straight lines and not along vertical ones, 
f the mode on plane A. as when o¢ = 0. The slope of the straight lines depends 
he on the magnitude of the delay, If the reference point 
in the sliding mode reaches the straight line L, with |5 | < 6, then, analogously to the case when o¢ = 0, it 
moves along the boundary Lg to the z-axis and enters a state of equilibrium(z = 0, 6 = 0, o = a¢). 
9) As in the case o¢ = 0, complex cycles and limiting cycles with partially sliding motion cannot exist in 
the system. 

In this way, the system has remained mainly unstable for the case ao, #0 too, The blind zone narrows 
stem the region of stability in Gg, which is now enclosed in the space between two cylindrical surfaces M and N, two 
om vertical planes C and D (Fig. 5), and above and below is bounded by the surfaces, the equations of which, in 

parametric form, are 
i 
21) 


= for |8|> 


Apart from this, the region of stability widens at the expense of the space from which the reference point 
arrives in the sliding mode regions of planes A and B with |5 | < max (6, 5g). Oscillations cannot arise in the 
system, as in the case when o¢ = 0. 


Comparing the system being studied with the system analyzed in [2, 3], we see that the plane's own damp- 
ing changes the structure of the phase space to such an extent that the system becomes mainly stable. Artificial 
damping cannot render the system mainly stable. It merely increases the region of stability by a small amount. 


5. Stabilization of Plane's Course by an Automatic Pilot with a Relay Servomotor 
Having a Blind Zone and a Loop in Its Static Characteristic 


We shall examine the stabilization of the course of a neutral plane, having negligible damping of its own, 
by means of an automatic pilot with a relay servomotor, having a blind zone and a loop in its static characteristic, 
The equation of the nonlinear function f(<) is, in this case, written as in (3). 


Phase space (z, 5 , o) consists of three parts (Fig. 8) which partially overlap each other. Part I (f(0) = 1) 
is bounded from below by plane A", the equation of which has the form 


A 


II ( f = 0) fs enclosed by planes A‘ and B* defined by the equations o=+(at 
(f(o) = —1) les below plane B*, defined by the equation 


A 


, Unstable limiting 


_ cycle 


_ 


Fig. 8. Phase space with og # O and A #0. 


The phase coordinates change in time in accordance with (4). Trajectories in parts I and III follow the 
vertical cylindrical surfaces which form the parabolas, defined in (5), at their intersections with horizontal planes, 
Projections of the phase trajectories on the planes z = Zp are determined from Eq. (6). Part Il fs filled with parts 
of phase trajectories which are situated in the planes 5, = 59 = const. Their equations have the form of (18). The 
joining of trajectories having the same direction takes place on the planes A and B' when & > 0 and on the planes 
B" and A‘ when 5 < 0. The sliding mode in the system is reflected in the motion of the representative (refer- 
ence) point between planes A‘ and A” or B’ and B". 


The transformations S*, characterizing the transition of the reference point (z», 59) from plane A* (5 < 0) 


. into the point (z4,54) on plane A* (6 > 0), has the form 


Bo = (C, T) =—¢, 


T? 
3, = 3, (6, T)=—C+ T,, 


4=46 


where T, is the time of the indicated transition. The 

{cycle transformation E*, which characterizes the transition 
Stable limiting of the reference point (z3, 53) from the plane A* (5 > 0) 
Cy into the point (zg, 53) on plane B", is written in the form 


‘Unstable limiting 


= 8 (C, T,) 


(27) 
= (6, Ty, T:) = + 


where T; is the time required for the reference point 
to cover the distance between the planes A* and B'. 


Fig. 9. Projection of phase space on the plane o = 0 ‘For the existence of a limiting cycle in the sys- 
for the case og #0, A+ 0. tem it is necessary to have 


8, = — — Zo. (28) 

i, 
From the first condition we find that T; = 2¢. 

With A = 0, parabolas L, and L, were the geometric 
positions of points satisfying this equality. Here these 
parabolas break up into four lines; curve Lj, lying in 
the left half-plane of A‘, curve Lf, lying in the right 
half-plane of A", curve Ly, lying in the right half-plane 
of B', and curve Lg, lying in the left half-plane of B” 
(Fig. 9). These curves are defined by the equations 


| 

From the second equation in (28) we find T;, 
La, taking into account that Ty = 2, 


Fig. 10. Diagram of point transformation for the case — — (30) 


o, #0, A#0., 


| 

| 
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n, 
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From the third equation of (4), setting f(o) = 0 in it, we find Ty, the time necessary for the reference 
point to cross from plane A” to the plane B': 


NE 


The curves ¥4(¢) and Y,(¢), shown in Fig. 10, have two points of intersection, Cay and {a,, which deter- 
mine two limiting cycles. An analysts of these curves, with the results of the studies set forth in the two preced- 
ing sections taken into account, shows that the limiting cycle with the small amplitude {a, is stable while the 
one having the large amplitude Ca, is unstable, 1. e., that the system, as before, with certain values of initial 
inclination, goes into runaway oscillations. Solving Eq. (30) analytically, we arrive at the equation 


1 


The roots of this equation, for small values of A, can be calculated approximately as in (12): 


(33) 


Thus, for small values of A, the range of initial conditions for which the system tends to oscillation practi- 
cally coincides with the small region of stability for the case A = 0, The amplitude of oscillations along the 6 - 
coordinate is directly proportional to the magnitude of the loop in the servomotor's characteristic. With increase 
in A the roots {g, and Ca, converge and for some value of A become equal, In this case, both the limiting cycles 
merge into one semi-stable limiting cycle, while the system becomes unstable for any initial conditions. 


The amplitude of the plane's oscillations along the angle of course can be determined by means of simple 
‘transformation 


4 


(34) 


An analysis of Expression (34) shows that ~, is mainly determined by the second term. For large N, when 
the first terms becomes large, Eqs. (33), which were used in the determination of , lose their validity. In this 
case €g, and {g, must be determined from graphical construction. The amplitude of oscillations along angle 
decreases with increase in o¢ (however, this is accompanied by a decrease in the range of initial conditions for 
which the system is stable) and with increase in the coefficients 1/a and B. 


The period of oscillation is T = 2(T; + T,), where T, and T, are calculated when ¢€ = Cay. 


As was pointed out earlier on, there can exist in the system a slipping mode which {s reflected in phase 
space by the motion of a reference point between planes A‘ and A” when 5< 0, and between B’ and B™ when 
6 > 0, The regions of the sliding mode on A‘ and A* and on B', B” coincide with the corresponding regions on 
planes A and B for the case when o, # 0, and A = 0 (Fig. 7). The difference lies in the fact thatfor A # 0 the 
amplitude and frequency of oscillation assume limiting values depending on A. 


A study of the phase space shows that there can exist in the system neither complex limiting cycles nor 
limiting cycles composed of slipping and nonslipping modes, 
In this way, for A # 0, the system oscillates for all values of its parameters. The amplitude of oscillation 


is determined by Eqs. (34). However, the system does not tend to break into oscillation for all initial conditions. 
As in the case of A = 0, the system breaks into runaway oscillations for certain values of the initial conditions. 


) 


2 
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The range of initial inclinations for which the system tends to oscillate approximately coincides with the region 
of attraction of the equilibrium position when og # 0, A= 0, 
SUMMARY 


1. The neutral plane, having negligible damping of its own, controlled by an automatic pilot with a con- 
stant speed steering servomotor, as a dynamic system is, for the case A = 0, manly unstable. 


2. The range of initial conditions for which the system tends to the equilibrium position (A = 0) or to os- 
cillation (A # 0) increases with increase in the object's inertia, steady feedback factor, coefficient of reaction 
along the derivative of the controlled coordinate, and with decrease in the servomotor's switchover time. 


3. The blind zone and the loop in the servomotor's characteristic decrease the range of initial conditions 
for which the system tends to an equilibrium position (A = 0) or to oscillation (A # 0). 
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ELECTROMECHANICAL CALCULATING DEVICE 


V. V. Gorskil 


(Moscow) 


The operating principles and description of the author's 
proposed electromechanica! calculating device (ECD) are given. 
Some results are given of experiments conducted on the simplest 
type of ECD, which may be used to plot Mikhailov curves, in- 
verse amplitude-phase characteristics of single-loop systems, 
and limits of D-splitting as a function of the amplification fac- 
tor, to make functional transformations of speed feedback, and 
to find real roots of algebraic equations. 


In the analysis of the stability and quality of automatic control systems, and especially in the utilization 
of frequency methods to this end, there often arises the need to calculate complex expressions. Such calcula- 
tions may be carried out on electronic miniature-scale setups and on digital computers, as well as on simpler 
electronic and electromechanical units, with an accuracy acceptable in engineering calculations. 


For this one may make use of, with some modifications, electronic and electromechanical {sographs — 
instruments for calculating the roots of algebraic equations [1, 2]. The simplest among these, although not the 
fastest, are the electromechanical calculating devices. 


1. Operating Principle and Description of the Simplest Type of an ECD 


The principle of operation of the is based on the property of synchronous generators, whereby they 
develop a voltage, depending in magnitude on the frequency (speed of rotation) and the excitation current, as 
well as on the frequency-dependent properties of electrical reactive elements (capacitors and inductors) which 


are supplied with power by synchronous generators. 


Figure 1 shows the main single line circuit diagram of connections in the transformation section and one 
channel of the ECD. In the lower part of the diagram are shown three synchronous generators SG;, SG3, SGs, 
which are brought into motion by the dc motor DPT, having a wide range of control over the speed of rotation. 
An emf, proportional in magnitude to the frequency (speed of rotation), is applied to the stator winding of syn- 
chronous generator SG, which has fixed magnets. A voltage U, = aww, proportional to the first power of the frequency 
is tapped off the potentiometer R,, which is connected through a rectifier to the stator winding of SG,. 


A voltage proportional to the square of the frequency is obtained by means of a quadratic function conver- 
ter [4, 5], consisting of a capacitor C, and resistor Rg. In fact, if the capacitive reactance {s much larger than 
the resistance over the whole range of frequency variation, then the current flowing in the circuit is I » Ey /x¢,* 
= kgCyw*, where E; = kg(w) is the emf of the first generator SG,. 


The voltage drop across the resistance Rg is proportional to the square of the frequency: 


U, = TR, = = (1) 
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For voltages proportional to the cube of the frequency, he second synchronous generator 1s used as well. 
Its field coil is fed by the stator winding of SG, through a rectifier and capacitance Cy. Since in this circuit too 
the capacitive reactance is much larger than the resistance Ry of the field coil of SG3, the excitation current 

flowing through it is proportional to the square of the frequency. If the magnetic system of SG, is not saturated, 
then this current produces a magnetic flux ®sG, proportional to the square of the frequency: 


where Ly 1s the self-inductance of the field coil in SGg. 

With the generator rotating with the same frequency w, flux ®5G, induces an emf E, proportional, in mag- 
nitude, to the cube of the frequency. 

The voltage Us = agw* is taken off the potentiometer Rs which is connected, through a rectifier, to the 
stator winding of SG,. 

A voltage proportional to the fourth power of the frequency is obtained by means of the functional conver- 


ter Cg~Ry. If the reactance of capacitor C, is much greater than the resistance R,, then the voltage drop across 
potentiometer R, will be proportional to the fourth power of the frequency: 


U, (3) 


Subsequent transformations are carried out by the functional converter Cs-R, , and by the third synchronous 
generator SG3, the field coil of which is fed through a rectifier and a capacitor C, by the stator winding of SG. 
Voltages, proportional to the fifth and sixth powers of w, are taken off R, and Rg: Us = agw®, Us = agu®. 

Further transformations may be achieved by the addition of more synchronous generators or by making use 
of known functional converters incorporating diodes, semiconductors, etc. In this way we obtain a series of volt 
ages, each of which is proportional to an integral power of the frequency. If we add to this series the voltage 
Vo = 49, which is independent of frequency and which is taken off potentiometer Ro, then the resulting sum of 

the voltages is equal (Fig. 1) to 


Uy = + + + + aot + + (4) 


It should be noted that similar transformations can be achieved using synchronous generators, without the 
use of intermediate functional converters, However, in this case, resultant voltages proportional to the first three 
powers of the frequency only can be obtained using the three generators. 


Equation (4) characterizes a whole range of problems which can be solved on such an ECD. 


Let us examine the accuracy of the conversion and the factors influencing it. The accuracy of the trans- 
formation depends on the relationship between the magnitude of the capacitive reactances and the resistances 
in the functional converters, on the voltage drops in the stator windings of the generators, on the magnitude of 
the armature reaction, on the saturation and residual magnetism in the magnetic circuits of the generators, and 
on temperature. 


The accuracy of the voltage conversion by means of a functional converter will be not less than 2%, if 

the capacitive reactance at the maximum frequency is five times greater than the resistance, To eliminate the 
effect of the voltage drop in the generator's stator windings, and the armature reactions, it is necessary to raise 
their power. At the same time the load current must not exceed the nominal current of the machine by more 
than 10-20%. The generators must work without saturation in their magnetic circuits, while the residual mag- 
netism must be eliminated every time before the commencement of work. Variations in the ambient tempera- 
ture practically have no effect on the accuracy of the transformations since in all the circuits the main imped- 
ance is the capacitive one. It should be noted that the apparent power of the synchronous generator is determined 
by its maximum voltage E and the permissible value of the excitation current in the following generator. The 
currents flowing in the load potentiometers and the functional converters (apart from those which are set up in 
the excitation circuits), inay be made as small as is desired. 


2. Experimental Assembly 


Figure 2 shows the main electrical circuit diagram of the transformation section and one channel of an 
ECD* model. The synchronous generators used were: SGy— an aeronautical speed-voltage generator type 2UG1- 
48 (two-pole, three-phase synchronous generator with fixed magnets, power 30 volt-amperes, voltage 75 v at 
3000 rpm), SGg and SGs— three-phase synchronous generators type FK 8,6-8 (eight-pole, 65 volt-amperes, cur- 
rent 1.05 amperes, voltage 36 v at 7500 rpm, frequency 500 cycles). 
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Since the voltages on the stator windings of SG, and SG, are relatively small, they are supplied with power 


through step-up transformers Tr, and Tr;, having one to five step-up ratios, so that the voltages across the resistances 
of the functional converters can be raised to 20-25 v. 


The values of the resistors and capacitors used in the circuit are shown in Fig. 2. Three-phase bridge recti- 
fiers made up of DG-Ts 24, DG-Ts25, and DG-Ts26 germanium diodes are used in the system, The generators 
are brought into motion by the dc motor DPT, which is controlled through a bridge rectifier by a laboratory auto- 
transformer type LATR-2. The same LATR-2 is used to demagnetize the magnetic circuits of SGg and SG,. The 
switchover circuit for demagnetization is not shown in Fig. 2. 


The constant voltage Vg» = ag is obtained from a regulated full-wave rectifier. In Fig. 3, the full lines re- 
present the experimentally determined relationships between the voltages and frequency, while the approxima- 
tions of these relationships are shown as dotted lines. From Fig. 3 it is apparent that the accuracy in reproducing 
the power functions aq, .. ., agw® on the ECD is quite acceptable for engineering calculations and designs. 


3. Examples of Calculations** 


The plotting of a Mikhailov curve and of the limits of D-splitting with amplification factor, the calcula~ 
tion of a transcendental function and the real roots of an algebraic equation, serve as illustrations, All the cal- 
culations may be carried out eithe in the true dimensional scale of the variable, or by going over to relative 
values of the variable and the corresponding transformed equation, The order in which the equation's coefficients 
are set up and the process of calculation are contained in what follows, Some value of frequency is chosen as 


* The model of the ECD was built by N. I. Karmatskii. 
** All calculations and drawings were done by R. F. Shepenina and E, A. Borisova. 
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a provisory unit, The magnetic circuits of SG, and SG, 
are demagnetized, motor DPT is started and the required 
frequency is established. The coefficients of the equa- 
tion are set up on potentiometers R,-R, with the help 
of a VTVM. DPT is stopped and SG, and SGy are de- 
magnetized, Then all the output terminals of the poten- 
tlometers are connected in series and are connected to 
the voltmeter. The signs of the coefficients are taken 
into account when connecting up the circuit. Readings 
are taken from the voltmeter and frequency meter as 
A the motor's rpm is slowly increased. For visual observa- 
oh tion, and to obtain a record of the output voltage as a 
ee function of frequency, an oscillograph with a memory 
3 tube may be used. 


Plotting a Mikhailov curve. Figure 4 shows the 
Fig. 5. variation of the real x(w) and the imaginary y(w) com- 
ponents of a Mikhailov curve as a function of frequency, 
for the equation 


p* + 2.5 p® + 3.32 p* + 3,725 p® + 1.67 p® + 0.389 p + 0.025 = 0. 


The calculated curves are shown dotted while the solid lines give the experimentally determined results. 


Plotting the limits of D-splitting with k + 1. Figure 5 shows the graph of the limit of D-splitting with re- 
spect to the parameter k + 1 for the equation 


— p> — pt — — 2.5 — 3p. 


The calculated curve is dotted while the experimental results are given by the solid curve. The range of 
variation of the relative frequency is shown on the curve. 
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Fig. 6, 


Calculation of e*. Figure 6 shows a positive exponential curve (dotted line) and the corresponding experi- 
mentally determined curve. The following series approximation was used here; 


at 


Extraction of real roots of an algebraic equation, Figure 7 illustrates the results of the extraction of the 
real roots of the algebraic equation 


y = — 1.2 + 0.68 — 0.08 —0.924 0.5 


by graphical means. The points of intersection of the curve y(x) with the abscissa correspond to the real roots 
of this equation. The dotted curve is the calculated one while the solid curve is the experimental one. 
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As was pointed out earlier on, the ECD or its separate circuits can be applied to the functional transforma- 
tion of speed feedback, the solution of complicated functions, etc, The range of problems soluble by the ECD 
can be considerably increased by introducing various nonlinear functional converters (combinations of capacitors 
and inductors). 
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4. Complete Schematic of the ECD 


The plotting of explicit frequency characteristics, curves of D-splitting interms of two parameters, and 
amplitude-phase characteristics is tied up with the solution of functions consisting of the ratio of two polynomials 
of the form 

+: ap + do 


Upon substituting p = jw, Expression (5) takes on the form 


(w) Ps (w) 
A (1) = (6) 


For the solution of (6) two channels are necessary, one of which will represent the numerator, the other the 
denominator, Moreover, the forming of the voltages in the second channel (denominator) is done by the same 
machines which are used for the first one. The complete structural schematic of the ECD is shown in Fig. 8. 
The voltages from the synchronous generators and the functional converters (not shown in the schematic) are im- 
pressed on potentiometers Ry-R, and Rj-Rj. The upper potentiometers are designated for setting up the numerator, 
the lower the denominator. 


Gi 

+ 


Fig. 8. 1) Numerator phase, 2) denominator phase, 3) squaring, 4) dividing 
circuit, 5) square root extractor, 6) frequency, 7) answer. 


The components P,(w) and Q;(w) are obtained by joining in one circuit all the voltages proportional to 
the even powers of w as well as the const2nt voltage, while to obtain the imaginary components P,(w) and Q;(w) 
the voltages proportional to the odd powers of w are joined together. Then these voltages are fed to the squaring 
circuit, are added in pairs and fed to the dividing setup. The extraction of the square root may be omitted if the 
instrument is graduated in the units corresponding to the square root of the quantity being measured. 


If an oscilloscope or an ordinary voltmeter is used as the measuring instrument then the extraction of the 
square root cannot be omitted. Ratio meters are used to measure the phases of the numerator and denominator. 
Division {s carried out by means of a ratio meter or any other arrangement. 


In this way the ECD, with two channels for setting up voltages, together with the squaring circuits and 
dividers, enables one to calculate amplitude~phase characteristics, plot frequency characteristics and curves of 
D-splitting in terms of one or two parameters, and, in general, to calculate the values of rational-fraction func- 
tions, 


SUMMARY 
1. The proposed electromechanical calculating device is sufficiently simple in construction, 


2. The preliminary experimental operations carried out on one channel of the confirmed the possi- 
bility of constructing a computer on this principle. 


3. The cited examples show that the computation accuracy {s acceptable for engineering calculations. 
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CIRCULAR TRANSMISSION OF REMOTE CONTROL SIGNALS WITH 
A COMBINATORIAL METHOD or\ CHOICE 


K. P. Kurdiukov 


(Moscow) 


Various methods of choice are investigated for time-separ- 
ated signals, and the possibility of circular transmission of signals 
using a combinatorial method of selection is demonstrated. Re- 
commendations are given for the use of interdependent pulse cod- 
ing together with a “stepping-switch" method of choice, and the 
effectiveness of such a code is determined. 


For the remote control of objects focused on a small number of massive points to be controlled (electric 
power stations, traction and step-down substations, etc.) by tele-controlling and tele-signalling (TC-TS) devices, 
the requirement is usually posed that there be guaranteed circular transmission of information signals concerning 
the objects controlled. Significantly less frequently is the requirement posed for circular transmission of the 
command signals controlling the objects, and special measures are even taken to guard against the execution of 
command signals containing more than one order. 


The great advantage of circular transmission of information signals is shown with particular sharpness in 
cases where it is necessary to re-establish quickly the picture of the controlled objects’ state after a disruption 
of the TC-TS device's functioning or during an emergency, and also for the control of connected or assembled 
objects for which a change of state of one induces a change of state in many others. Thus, for example, a devia- 
tion of the supplying lead-ins to the traction substations of the Moscow subway engenders up to 39 different switch- 
overs, necessitating the same number of tele-signals from the dispatcher's panel [1]. 


In these conditions, circular systems, to a large degree, meet the requirements of maximum clarity of se- 
presentation, on the dispatcher's panel, of the occurrences at the controlled points, and minimum duration of 
signal transmission for switch-overs (not isolated ones) requiring the attention or intervention of the dispatcher. 


Precisely this requirement for circular transmission of information signals was predetermined in a host of 
cases making use of the distributive method of selection in TC-TS devices with time-separated signals, 


Ordinarily, for remote control connections in remote-controlled industrial and transport enterprises, inde- 
pendent channels are used (in the remote control circuits) and, therefore, the question of transmission effective- 
ness and rapidity of action of TC-TS devices was not sharply raised until quite recently. However, with the in- 
crease in speed of the controlled and regulated objects themselves, and with the increase in the number of objects 
included in the channel, this question acquired ever greater importance and forced a reconsideration of the pre- 
viously accepted solutions. 


With respect to speed and effectiveness of transmission of individual signals, the distributive method of 
selection is inferior to other methods, since it requires the expenditure of one code element for each signal 
which could be transmitted. In other words, with the distributive method of selection, the code must contain 
a number of elements (pulses or intervals) equal to the number of signals, although, as stated above, in any one 
cycle of tele-control, only one error may be transmitted. 


TABLE 1 The formula for the code for a distributive method 


of selection has the form: S =n, where § is the total 
nies unitate number of signals and n {s the total number of code ele- 
Signal {Code elements ments. 
no. 1! 2|3]4 | 5] 6 The circular possibilities of a code may be esti- 
mated by the coefficient of circularity, 5 = s/S, where 
1 1/0/)0/0/0/0 © s is the number of signals which can be transmitted in 
one cycle, For the distributive method of selection, 
4 6 =1, 
5 
6 0/0}0\0)0}14 The distributive method, in its pure form, is rarely 
used (for example, in the TC-TS device type VRT~-53). ' 
In order to increase transmission effectiveness, the distri-~ 
TABLE 2 buted group method is ordinarily employed where, for 
example, the first selecting pulse of the code chooses 
Distributed Group Method the group of objects and the second code pulse chooses 
the object in this group. Thanks to this, the code may 
Grou Object contain fewer elements than there are signals. The for- 
Group | Signal} selection | selection mula for such codes has the form: S = ny, where ny 
no. eel Er is the number of elements for p selection, equal to 
Goce sine the number of groups, and ng hore number of elements 
2; | for object selection in each group. 
‘ 110 1/010 For such a code, the circularity coefficient is 
2 
s 
4 
It is well known [2] that {s obtained with 
the condition, ny = ng= vS, 
7 |0/0/;14/)/1}0/0 The use of the distributed group method of selec- 
tion is sometimes brought about by the presence of some 
executive points which,in this case, are considered as ) 
independent groups of objects. The first degree of selec- . 
TABLE 3 tion is then the choice of such isolated groups. 
Combinatorial Method of Selection from One Set It is necessary to bear in mind that the introduc- 
tion of each subsidiary group and degree of choice lowers 
Signal _Code no, the circular possibilities since, in circular transmission, 
no. i 2 3 | «| 5/6 only signals from one group may take part. Thus, for 
example, with two groups, circular transmission can en- 
1 4 1 0 |0/;0/0 compass only half of all objects, In this case, the coef- 
ficient of circularity 6 will equal one half. | 
~ For the analysis of methods of selection, it is con- 
6 0 1 1 |ololo venient to represent the signals in the form of binary 
7} 0 i 0 }1/0)0 numbers, and the set of signals in the form of tables. 
Here, the selecing code element will correspond tol , 4 
10 0 0 4 |414;0)0 and the motion element to 0, 
11 0 0 1 |0)1/0 
12 0 0 1 |0;0)14 In Tables 1 and 2 are given the sets of signals for 
the distributive and the distributed group method of 
15 0 0 0 |O;4];1 selection, constructed by decreasing binary numbers for 


n =6, 
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TABLE 3a 


Circular Combinatorial Method of Selection from One 
Set ("Stepping-Switch” Method of Selection) 


Group | Signal Code no, 
no. | no 
1 | 2 
1 i 1 0;0;0;0 
2 1 0 1;0/;0/0 
i 3 1 0 0;1;0/0 
4 1 0 0/;0;1/0 
5 1 0 
6 0 1 1;0/;0/0 
2 7 0 1 O;1;0;)0 
8 0 1 0;0;1/0 
9 0 1 
10 0 0 1 0; 0 
3 11 0 0 1;0/;1/0 
12 0 1 
4 13 0 0 O;1/;1)0 
14 0 0 
5 | 5 | o | o jololala 


Fig. 1. Circular code for the “stepping-switch” me- 
thod of selection: a) choice of the second, fourth, 
sixth and eighth objects from the first group, b) choice 
of second and fourth objects of the sixth group, c) the 
selection of the one object in the tenth group. 


In Table 3 is given the set of signals for the com- 
binatorial method of selection from one set with a choice 
of two, also for n = 6. The formula for a code when this 
method is used has the form: §S = —- = C2, where m is 
the number of selecting elements of the code. 


Here, each object is selected by means of two 
code elements, exactly as in the distributed group me- 
thod, but there is no longer a constant division of the 
code elements according to the level of selection. How- 
ever, in analogy to the distributed group method, it may 
be considered that the first selecting element selects 
the group of the object, and the second selects the ob- 
ject from its group. Using this idea, one may divide 
the selecting code elements into groups, as shown in 
Table 3a. 


Considering all the combinations of Tables 3 and 
3a from this point of view, one may easily convince one- 
self that, within the second level of choice, i. e., the 
level of object selection, it is completely possible to 
carry out circular transmission of signals without vio- 
lating the principles of combinatorial selection. The 
boundaries of the selection groups are shown by the hori- 
zontal lines on Table 3a. 


Due to the variable capacity of the groups, the 
circularity coefficient is different for the various groups, 
and varies from 6 =(n-1)/S tod = 1/S for noncircular 
transmission, characterizing the last group, which con- 
tains only one object. 


This is also clarified by Fig. 1, on which, as exam- 
ples, are given the codes for two-level selection by the 
“steppoing-switch" method, where the selecting marks 
are agreed to be taken by increasing pulse amplitude. 


It should be mentioned that, for the same number 
of code elements, the distributed group method cannot 
guarantee as high a circularity coefficient as the most 
capacious group using the combinatorial method since, 
even with just two groups, the circularity coefficient 
5 , for the distributed group method, equals (n — 2)/S. 


Using the combinatorial method of selection from 
one set with a three-fold choice (analogous to a three- 
level application of the distributed group method), one 
may consider that the first selecting code element selects 


the group of the object, (for example, the group of isolated objects, i. e., the controlling points), the second ele- 
ment selects the subgroup of this group, and the third element selects the object in the subgroup. 


Extending this to arbitrary combination of n things taken m at a time, one may write the general expres- 


sions for the number of signals, the number of groups and the number of objects in the groups, and also the cir- 
cularity coefficient for a given number of code elements for such a circular combinatorial method of selection 
from one set which, due to the way in which it was developed, was called the “stepping-switch” selection method.* 


* The “stepping-switch” selection method was laid down, and developed, by K. P. Kordiukov and L, I. Kazan- 


skii. 
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The total number of signals (the code formula) is S = Ch’. The number of groups is g = C™5", The number of 
signals, s, in the group varies from n—™m+1tol, The circularity coefficient, 5 , varies from 


n—m-+i 1 


Figure 2 gives curves characterizing the selection methods considered here, The total number of signals 
is laid out along the abscissa and along the ordinate axis, the number of pulses which must be contained in the 
code in order to transmit these signals. This quantity can also characterize the duration of transmission, if the 
duration of the code elements is known and if the difference in duration between the selection and motion ele- 
ments is ignored. The dotted lines characterize the number of signals which may be transmitted circularly in 
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one cycle where, for the “stepping-switch” method, the number of signals of the first group is stated. On Fig. 2, 
curves 1 and 1° (shown as separate although they actually coincide) are for the distributive method of selection, 
curves 2 and 2° are for the distributed group method for two groups, curves 3 and 3° are for the distributed group 
method with the number of groups, g equal to vS, 4 and 4° are for the combinatorial method of choice by two 
elements, and 5 and 5° are for the combinatorial method with a choice by m, where 


n-+-1 


m= or m= —— 
2 


3 


Figure 3 gives curves for the dependence, on the total number of signals, of the number of counting, group- 


ing and composing relays necessary in the decoders for the various selection methods. This number fs calculated 
from the formula, K = n + g. 


ab 


az 


Fig. 4. 


Here, curve 1 is for the distributive method, curve 2 for the distributed group method for two groups, curve 3 
for the distributed group method with the number of groups equal to g = VS, curve 4 is for the “stepping-switch" 
method with a selection by m and curve 5 for the “stepping-switch” method with a choice by two. 


The number of executive relays for this was not counted, since the number of them does not depend on the 
method of selection, and is always proportional to the number of objects. Likewise, the number of protective 
relays, the number of which depends on the structure and design of the device, was not calculated. 


TABLE 4 Figure 4 gives the curves for the dependence of 
Distributed Group Method the relative number of selecting and distributing relays 
per signal on the total number of signals. The nomen- 


lature is the same in Fig. 4 as in Fig. 3. 
Group | _ 1st degree | 2nd degree — ig ig 
no. ag Code elements It is clear from Figs. 2-4 that the "stepping-switch" 
a Pa n, Pr | pa | method of selection, somewhat inferior to the distributive 
ttialololololilolo method with respect to circularity, requires the lowest 
4 number of relays and decoders. 
1ol1lololilololo It is necessary to state that, in using the “stepping- 
23 410/1/0/0/0/0/11/0 switch" method of selection, it is necessary to take cer- 
tain additional measures to prevent the execution of 
3 4 false commands, arising from defective devices or from 
noise on the line. If, for example, the first selecting 
h 7/0/0/0/1)1/0/0)0 code element, due to noise, was lost or somehow in- 
completely reproduced then, using circular transmission 


of the signals, i. e., with some object-selecting element 


aps 
it 
| 
6 8 20 90 100 209 $ 


TABLE 5 in the code, the first element would seem capable of 


é » lecting an incorrect group, and the second element 
Method 
Sine; ae amines of selecting the incorrect object from this group. 
Group /Signal) _Code elements For protection against this eventuality, there might 
Ps |ps ns | Ps | m | | be recommended the use of different type pulses or code 
oP elements according to the level of selection. Thus, for 
4 1 example, with a time-pulse code, for the first level of 
3 11/0}0);0)/0)14];0]0 choice the length of the pulse might be used and, for 
4 |1/0/0)}0)0)0)0)14 the second level, the length of the pause, or vice versa. 
In this case, the loss of the first selecting sign would 
5 |0/1/1,0)0/0/0/0 lead to an illegal code, i. e., to the prevention of fail- 
2 6 0} 4] 0] 0] 0 
7 
It is possible to make an even fuller use of the 
selecting signs when part of the group is chosen by one 
10 lololailolololol4 sign and the other part by a different sign. The corre- 
sponding object is also chosen by the two signs but in 
4 414 |O;0}0/1/1])/0])0)]0 a completely contrary way to that in which the given 
12 |0/0/0)1/0)0)1)0 group was chosen. This might be called the crossover 
‘ 13 lololololalalolo method of using two selecting signs. The signs themselves 
144 |O|O|O;}O}4/}/0]10]4 might be interdependent, i, e,, be such that, under the 
7 action of noise, can be converted one to the other as 
| 0 eyo occurs, for example, in the time~-pulse code using the 
7 | 16 JOjJojo}o|ojo|4)4 length of pulse and pause. The cross over method of 


using them guarantees complete protection from the 
execution of signals thus distorted. 


It was found [3] that the total number of signals (i. e., the code formula) and the number of groups for the 
crossover use of pulses and pauses are expressed in the following way: 


for the distributed group method 
5Sp.n= Pia + MP2, = Pr tM; 
for the "stepping-switch" method with selection by two 
+Crand g =Cy+ 


for the “stepping-switch" method with selection by m (m an even number) 
Soe 
g= 


er a+ 


for the "stepping-switch” method with selection by m (m an odd number) 


_ 


where p is the number of pauses, n is the number of pulses and m is an odd number. 


Tables 4 and 5 give the set of signals for the crossover method of using pulse signs for, respectively, the 
distributed group and “stepping-switch” methods of selection, for p = n = 4. 


The curve, S= f(p, m), for the crossover method of using pulse signs with a “stepping-switch” method of 
selection is shown in Fig. 2 (curve 6). It should be mentioned that, with the use of the crossover method and the 
“stepping-switch® selection method, the number of signals grows to more than twice (as occurs, for example, for 
the distributed group method) in comparison to the separate use of one or another sign. 


SUMMARY 


The “stepping-switch” method of selection may be considered as a variation of the distributed group method 
in which there {s no constant division of the code elements by level of selection (apart from the m — 1 first and 
successive elements), or as a variation of the combinatorial method of selection from one set in which is guaran- 
teed circular transmission of signals in the group selected by the first m— 1 selecting code elements. 


The “stepping-switch” method of selection combines in itself the positive qualities of the distributive me- 
thod (circularity) and the combinatorial method from one set (higher transmission effectiveness). 


The “stepping-switch” method of selection by two requires the least outlay of relays and decoders, 

The “stepping-switch” method of selection together with the crossover method allows the use of two inter- 
dependent pulse signs, which, with full protection of the signals, more than doubles the effectiveness of transmis- 
sion. 
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CONCERNING ONE TYPE OF COMMUTATION CIRCUIT 


V. G. Lazarev and Iu. L. Sagalovich 


(Moscow) 


One type of commutation circuit is investigated — a cir- 
cuit for converting binary numbers for every possible permuta- 
tion of their terms. The method of synthesis, based on the use 
of the theory of permutation groups, provides the possibility of 
obtaining optimal circuits for a small expenditure of trouble on 
their components. The formula will be provided for the com- 
putation of the number of contacts necessary for the construction 
of the commutation circuit considered. 


Together with the general methods of synthesizing relay-contact circuits [1-4] there exist particular me- 
thods for individual types of circuits, allowing significant simplifications and increases in speed of the process 
of planning these types of circuits. The present notes on the construction of a commutator present one such 
method, which makes use of certain results from the theory of groups.* 


The problem of constructing a commutator arose in connection with the development of a machine for 
the synthesis of contact circuits [6], modeling the method of cascades [3] (the graphical method of [4)). 


It is well known that the number of contacts in a circuit which has been synthesized by the method of 
cascades, and certain other of its parameters, depend on the order of synthesis [3], i. e., on the order in which 
the selected elements follow one another. If the number of such elements is n then, to obtain the optimal (for 
the given method of synthesis) circuit, it is necessary , in the most general case, to examine all n! circuits ob- 


tained as a result of carrying out all possible permutations of the numbering of the variables in the decomposi- 
tion into constituent units of the function to be realized. 


The equivalence of the problem of permuting the numbering of the variables with the problem of substi- 
tuting binary numbers becomes obvious, amounting only to the decomposition of the function to its constituent 
units whereby each nonnegated variable is treated as a 1, and each negated variable as a 0. 


We remark that an arbitrary permutation of the terms of an n-place binary number does not change the 
number of ones and zeroes in the number. Therefore, all 2" n-place binary numbers may be divided into n + 1 
classes and, consequently, all commutators may be divided into n + 1 sections (Fig. 1), adopting the number r 
for that sector assigned to binary numbers containing r ones, The class numbered r contains Ch, binary numbers. 

For an arbitrary nonidentical permutation of terms, the 


44 binary numbers, in some fashion, change places within 
ao 4 ah ee Li bee the classes, since there is no pair of terms with the same 
WO\ Wel Nor elements in all the numbers, 
IT Thus, to each permutation of order n of the terms 
/ n Af, Ar / there corresponds one and only one permutation, of 
Fig. 1. order Ci, of the binary numbers themselves. 


* A particular case of the circuit presented was studied in (5). 


The following assertions hold. 


1) The permutations of the binary numbers, obtained by the permutation of their terms, form a group iso- 
morphic to the symmetric group Sn [7]. 


In fact, let the binary number Aj, under the identity permutation of terms, with the ordinal numbering 
1, 2,...,m, have ones in the binary places; 


* Hig, Bir, ion 1, 2, Cor 1< aj qn. (*) 
The first subscript shows that the given “one” is to be found in binary number A;, while the second sub- 


script gives the ordinal numbering only for those binary places in the given number which contain “ones" The 
collection (* ) is the set of all possible combinations of the n digits, 1, 2,...,m, takenr at a time and, there- 


fore, any symbol oj; is repeated in exactly C'”" sets, 1. e., the subscript 1 assumes Cl"! values simultaneously. 


n-1 


We consider the permutation of terms 
( 1 
= 
. 


(the first subscript has been omitted). It says that the “one” which was at position p now goes over to position 
Op» i, e., that the “one” which was found in one of the os binary numbers is now, generally speaking, to be 


P 
found in another of the Cr”! binary numbers. In this case, by repeating each 1.) cri times, and rear- 


ranging them in the corresponding way (so that the permutation itself is not changed), we can give the permu- 
tation s in the form: 


Ch 


Permutation s; implies the permutation 


In exactly the same fashion, the permutation 


n n 


implies the permutation of the binary numbers 
A, Ai, - A, 


Here, the product of the permutations Sj and Sj will be the permutation 


‘ 


Ace 
n 
\ 


since precisely this is implied by the permutation 


$ = $s; = > 


n n 


Thus, between the permutations s of the terms and the permutations $ of the binary numbers is established 
a one-to-one correspondence and, if permutation sj corresponds to permutation S; and permutation 8; Corresponds 


to permutation S;, then the permutation ss; corresponds to the permutation S{Sj, which is required for there to be 
an isomorphism. 


Thanks to the existence of this isomorphism, the 
w set of all actual permutations of order Cf, of the binary 
numbers constitutes a group H and 


| 2) Every transposition (ij) < Sn corresponds to 

a permutation H, which decomposes into a 

t product of independent binary cycles. This assertion 

Fig. 2 follows immediately from the fact that the correspond- 
= ing elements of the isomorphic groups are always of 

one order. 


It is easy to compute the number of independent binary cycles of the permutation S;,. 


Indeed, under the transposition (ij), the terms are changed only in those binary numbers which contain either 
0 and 1 or 1 and 0 in the transposed terms. The one and the other will be possible only the number of times that 
r— 1 “ones” can be placed in the remaining n— 2 binary places, i. e., Cr_4. This is the desired number. The 


converse is also obvious; all products of Ch-4 independent binary cycles of the group H correspond to a trans- 
position of the group S,. 


3) It is well known that the set of 4 transpositions is a system forming the symmetric group ©, [7]. We 
conclude from what was stated sbove that the corresponding products of Cl~4 independent cycles is also a sys- 
tem forming a group H. 


After having made these remarks, we turn to the construction of the commutator or, more accurately, to 
that portion of it which permutes a n-place binary number containing r “ones.” We shall call this the r-th sec- 
tion of the commutator. The r-th section of the commutator takes the form of an orientated (Ci, Ch)-pole 
(Fig. 1), the input and output of which represent the binary number A; (i= 1, 2, .. ., Cj). Conduction between 
input A; and output A; must occur each time that the permutation modeled by the commutater transforms the 
number A; into the number A; , the permutation itself being of the terms of the numbers. 


By a “cross section" of the commutator we shall mean those lines passing through all contacts of one and 
the same relay modeling transposition (ij). Due to the noncommutativity of multiplication of permutations, the 
cross sections must not intersect. Since a relay isa two-position device, the contacts of this relay, by virtue of 
2) above, must be connected in independent blocks, as shown in Fig. 2; for each cross section there are, conse~ 
quently, such blocks, 


It is easy to show that, with any arbitrary ordering of the C? transpositions, with their permissible order of 
multiplication maintained, it is possible to obtain the entire symmetric group* 6G,,, and, by virtue of 3) above, 
also the group H. 


* We mention that, in [5], it was only proven that all C? transpositions are sufficient for obtaining group G,, . 
but it follows from the method of proof that all C2_, transpositions which do not vary symbol n precede all the 
transpositions which do vary that symbol, 


This means that the cross section of the r-th commutator section may be laid out in any order. Using the 
results of [5], one can reduce the size of the cross section of the r-th section (and, consequently, of the entire 
commutator). 


It is easy to see that, for each binary number, there is a double for it, in the sense that wherever there is 
a “one” in one number there fs a “zero” in the other. Therefore, the r-th and the (n— r)-th section of the com- 
mutator are implemented identically, and are isolated from each other. 


If n is even, the (n/2)-th section cannot be decomposed into two isolated parts, since the double of the 
binary number with the same quantity of “zeroes” and “ones” may be interchanged with the number itself. 


The number of switching contacts of the entire commutator is given by the formula 


N = 2m C33, 


r=1 


where m is the number of cross sections. 
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THE APPLICATION OF A THERMOELECTRIC PROPORTIONAL-~PLUS-INTEGRAL 
CORRECTIVE DEVICE FOR IMPROVEMENT OF TWO-POSITION 
TEMPERATURE CONTROL 


A. A. Kampe-Nemm 


(Leningrad) 


An automatic control method which is based on applica- 
tion of a thermoelectric corrective device and combines the pro- 
perties of two-position and proportional-plus-integral control is 
treated. The construction of the corrective device is described. 
Experimental results show that the proposed device permits a sub- 
stantial improvement of the two-position control process. 


The feedback-type corrective devices applicable in relay systems of automatic control are described in 
[1-13]. As has been shown in these works, this kind of corrective device can substantially improve the control 
quality of the relay systems with constant speed servomotors, as well as of the two-position control systems. The 
feedback systems in the form of thermoelectric corrective devices and their ap- 
plication with two-position temperature controllers were described in particular 
in [4-6], It has been shown that this kind of device allows a sharp decrease in 
oscillation amplitude of the controlled variable, but at the same time they make 
the system static with a residual unbalance, 


As an addition to the articles mentioned, there follows below a descrip- 
tion of a device* which makes it practically possible to eliminate the residual 
unbalance and gives a system of the proportional-plus-integral type. 

In the diagram of Fig. 1 the following nomenclature is used; 1 is the ob~ 
ject of control (an installation with an electric heater), 2 is a two-position-type 
control unit, 3 is the corrective device, H is the heater (of the object), T is a 
thermocouple, K is an electric contact. 


Conductor C is the basic element of the corrective device; it is made of 
a material different from that of the conductors connecting thermocouple T with 
unit 2, The ends of the conductors T, and T, are heated by heaters H; and Hy 
which are turned on and off with heater H. Heating combination H,~T, has a 
faster response, while H,-T, has a slower response. When contact K ts closed 
thermal emf's, which add algebraically to the thermal emf of thermocouple T, 
are developed at the ends of the conductors T; and T,. 


Junction T, serves to create a pulsed mode of operation, while junction T, 
serves to eliminate the residual unbalance, Due to slower response of junction 
Tz it undergoes smaller temperature changes than does junction Ty. But since 
the heaters of both junctions are turned on and off for identical intervals of time, 


* Authorship Certificate No. 104760 issued to A. A. Kampe-Nemm and A, N. Trusov. 
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the average temperatures of the junction become the same after the transients die out; therefore, the average 
values of emf developed by them are cancelled out by each other and the residual unbalance is eliminated. 


For the corrective device to operate property, the extreme values of temperature (emf) to which both junc- 
tions could be heated if the heaters were energized for a longer period of time, should be the same. It is also 
essential that the time constants of both heating devices be sufficiently different from each other. 
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Fig. 3. 


The described temperature control circuit has been experimentally tried out on an electric muffle furnace. 
The construction diagram of the corrective device used is shown in Fig. 2. Element I, with a faster response, con- 
sisted of two mica plates 1 and with a junction of two conductors, copper and constantan, with T, squeezed be- 
tween them, Nichrome heater 5 was placed outside. Element II, with a slower response, differed from the faster 
element by the fact that an aluminum plate 6, shown by a broken line, was put between the mica plates. It was 
in contact with junction T, of two conductors, constantan 3 and copper 4. By means of binding posts A and B a 
voltage was applied to the heaters, while binding posts X, Y and Z served to connect the device with the thermo- 
metric circuit. In the faster element 1 a sheet of foil 7, shown by a broken line, has been placed between mica 
plates. It served to secure a more uniform heat transfer from the heater to the thermocouple. 


Figure 3 is a graph of automatic temperature recordings obtained during the tests. Prior to instant t, the 
control, with the corrective device taking part, was in effect. The furnace. doors were closed before t; and after 
t, and were open during the interval between t; and tg. As can be seen from the graph, the controller {s well 
able to take care of the control even under such strong disturbances. For comparison there is given, on the right 
of Fig. 3, a graph of the two-position control without the corrective device, where in the interval between t, and 
t, the furnace doors were closed and after t they were open. In both instances the control process is accornpanied 
by large oscillations. 
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SUMMARY 


1. The application of the described thermoelectric corrective device allows a considerable decrease in the 
oscillation amplitude of the controlled variable, as compared with the conventional two-position control. Under 
undisturbed conditions the control curve becomes practically a straight line. 


2. The advantage of the proportional-plus-integral corrective device over that of proportional-pulsed 
(static) type consists in its ability to eliminate the residual unbalance and to bring the controlled variable to 
the assigned value with sufficient accuracy even under strong disturbances. 


3. Proportional-plus- integral thermoelectric corrective devices are simple to build and can be used as 
additions to two~position controllers to improve their performance. 
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IMPROVING THE QUALITY OF TWO-POSITION CONTROL 


A. I. Cherepanov 
(Molotov) 


The operating principle and the construction of a correc- 
tive device proposed by the author for a two-position controller 
are considered. Some results of a graphic analysis of the process 
of two-position control with correction are given. 


As is known, the quality of two-position control can be substantially improved by means of special correc- 
tive devices, as for example, those described in [1, 2]. 


In this article it is shown that the precision of a two~position controller can be increased by introducing 
a unique kind of interruption into the control process by means of a special corrective device of the time relay 
type. In the case of an extremely unsymmetrical two-position control process this kind of correction is simply 
realized by means of an electromechanical device, 


| a which would serve to limit the amplitude of only one 
10 (for example, positive) deviation of the controlled para- 
‘ meter, 


Let us note that the unsymmetrical mode of opera- 
tion of two-position control takes place in many prac- 
tical cases, as for example, in relatively low tempera- 
ture control in thermostats, drying racks, muffle furnaces 
and others, 


The basic diagram of an electromechanical cor- 
rective device is shown in Fig. 1. 


The corrective device has a supplementary relay 
P with two pairs of contacts: P-2— normally closed, 
and P-1— normally open, Relay P is connected in 
parallel with heating element HE and is energized when 
control contact K is closed, The contacts of relay P 
send the controlling pulse to the time relay, which con- 
sists of an electric motor M driving two discs: a basic 
Fig. 1. disc D-1 and a supplementary one D-2. These discs 
control contacts K~-1 and K~2, respectively, Discs D-1 
and D-2 are affixed on the same axle and their cams are so constructed that adjustment of different timing inter- 
vals is possible. 


Figure 2 shows the process of two-position control for the case of a controller without a nonsensitive region 
(curves AS and A®%, . ) and a process in the same system with the corrective device operating (curves A 9 and 
4 9 ¢,.)- In Fig. 2 the following nomenclature is used: A is the deviation of the controlled parameter (tem- 
perature) of the object, AS, . is the temperature deviation measured by a sensitive element, A 9,4) is the 
amplitude of positive deviation, A.) is the amplitude of negative deviation, T,,, is the time on, T off is the 


time off, T is the period, t, is the instant the inflow {s shut off in an ordinary control, The same quantities de- 
signated by primes correspond to a system with the corrective device, Qin is the heat inflow, Q,. is the heat 
outflow, At is she delay time, At" fs the “residual” delay time with the correction process taking place, tj is the 
stopping time of the discs after the first revolution. 


The operation of the corrective device is based upon the independent shutting off of heat inflow to the ob- 
ject an assigned time interval after the controller contact is closed. 
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Let us consider the action of the device from the instant when the object temperature starts to fall. The 
temperature measured by the sensing element (A 9, ¢ ) lags the temperature of the object. Att = 0 the controller 
contack K closes and energizes the supplementary relay P, which in turn energizes the electric motor of the 
time relay by closing its contact P-1. At the same instant the object temperature begins to rise while the tem- 
perature of the sensing element continues to fall. After an assigned time tj contact K~1 is opened by the cam 
on disc D-1 and the energy inflow to the object is stopped. At this instant the object temperature begins to fall 
after having passed over the assigned value and having reached its maximum A Om+)- At instant ty the sensing 
element temperature reaches the assigned value and contact K closes, It is obvious that at thisinstant or some- 
what later contact K-1 can be closed again. It is closed by cam on disc D-1 at %. After ty the time relay con- 
tinues to operate, since relay P has blocked the motor through its contact P-2 and the closed contact K-2. At 
some instant t3 contact K~2 is opened by the cam on disc D-2 and the motor stops. The corrective device is 
ready for its next cycle. 


By utilizing the simplified method of graphic analysis for computation of two-position control processes 
presented in [3], the expression for basic parameters of the corrected process can be obtained. 


We are introducing these formulas in their final form. 
The amplitude of positive deviation is 


Qin — 2 


where C is the thermal capacity of the object. 


The amplitude of negative deviation, as can be seen from the graph, remains the same: 


The time on is 


Ton 


The time off is 


Q 
= — Ata (2) 


= (hs Au’ 
+). 


Qout 
Qout 


(3) 


(4) 


The period of oscillation of a system with the correction has the form 


= 
Al 
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mity, 


+ Sin 


(5) 


Denoting the inflow-to-outflow ratio by p = 


= — and introducing the concept of degree of cor- 
out 
rection k = rd , from Expression (5) we can obtain for 


At =1 sec 


n[k + (n—1)| 


Figure 3 shows the relation between the oscilla- 
tion period T’ of a system with a corrective device and 
the inflow-to-outflow ratio n for various values of the 
degree of correction k, The curves show that the oscil- 
lation frequency increases with increasing degree of 
correction (i. e., the precision of control), but for higher 
values of k the frequency increase slows down. 


Parameters of a corrective device under steady 
state conditions can be selected by means of the for- 
mulas derived above, In doing this, it is necessary to 
begin with the assigned precision of control, {. e., the 
admissible value of A%,,,). Obviously the degree of 
correction is determined as 


The time when contact K~1 is closed 1s obtained from Formula (3) and the expressions for n and k: 


k(n — 1) 


The time when the contact is open is found from the formula 


k— 


1 
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For the minimal time value of the new closing of the contact, we have 


At, 


n—i 


According to the curves of Fig. 2, the time of one revolution of the discs under steady state conditions must 
satisfy the inequality 


or 


n n [k + (n— 1)] 


Without going into the details of transient processes in a two-position control system with a corrective de- 
vice, we shall only note that in the case of prolonged constant disturbances on the outflow side, there can appear 
complex periodic motions in the system, In such a case, because of the forced temporary circuit break, there 
will be an excessive decrease in object temperature (the disruption of control) for values of inflow-to-outflow 
ratio somewhat greater than unity. Let us call this value nc. Judging from these considerations, the instant of 
the new closure tj should be close as possible to ty, Sometimes it is possible to choose t3 even somewhat less than 
ty. The time of one disc revolution, however, should be chosen to be as high as possible and little different from 

T*. When the parameters of the corrective device are 


The instant the correction is turned on selected, the nature of the disturbances of the system 
should be considered. 
w The expression for the critical value of ratio 
can be obtained from the time relation of the closed and 
$ opened conditions of contact K-1 in one complete disc 


revolution, As a result the formula* 


(n — 1)] 
cr n? + (k— 1) 


Fig. 4. 
is obtained where n is the nominal value of the inflow-to-outflow ratio (under steady state conditions). 


The experimental check on the proposed method was conducted with a system where the controlled object 
temperature was that of a drying cabinet. 


The temperature was controlled and recorded by means of an electronic potentiometer EPD-12 and a chro- 
mel-copel thermocouple without a jacket. 


The process of ordinary two-position control, the process of two-position control with correction and the 
parameters of the corrective device were calculated, 


The given design data are; the power rating of the cabinet N = 1 kw, the thermal capacity C = 4 kcal/deg, 
delay time At = 90 sec, inflow-to-outflow ratio n = 6. 


The following basic parameters of the two-position control process were obtained by calculation: positive 
deviation amplitude A 9,,(,) = 4.5°C, negative amplitude A9,,_) = 0.9°C, period T = 10.8 min. 


It was necessary to secure +1°C accuracy of the temperature regulation in the cabinet. 


The basic parameter values obtained for the process with correction were; positive deviation amplitude 
A Sin(+) = 1°C, degree of correction k = 4.5, negative amplitude A 9.) = 0.9°C, time on T,, = 0.64 min, 
period T* = 3.84 min. 


The accepted parameters of the corrective device were as follows; time off tj = 0.7 min, time on & = 
= 1.7 min, time of one revolution tj = 3.5 min. 


The experimental record of this example is shown in Fig. 4. 


* The formula is derived for maximum possible time of one revolution, 1. e., for ts = Ty. 
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LETTERS TO THE EDITOR 


ON THE RELATIONSHIPS BETWEEN THE ERROR COEFFICIENTS AND 
THE AMPLITUDE AND PHASE CHARACTERISTICS OF LINEAR 
REPRODUCING SYSTEMS 


Vv. G. Vasiltev 


(Moscow) 


Expositions of the connections between the error coefficients and the other parameters and characteristics 
which describe reproducing systems have appeared in a number of works [1-3]. In particular, the question of the 
interrelationship of the error coefficients and the system's frequency characteristics (imaginary, real and ampli- 
tude) was dealt with in the paper of A. M. Boev. 


However, in A. M. Boev's paper the interrelationship of the error coefficients and the system's phase charac- 
teristic remained unexplained. The present note establishes the character of this relationship. 


1. It is well known that the behavior of the system's transfer function along the imaginary axis of the p- 
plane is unambiguously determined by the system's amplitude and phase characteristics, M(w) and ¢(w). The 
following equations express the interrelationships of these characteristics with the system's transfer function, K (p): 


K (4 j@) = M (w) 


(la) 

M (a) = V K (ja) K (— je), (1b) 
K 

(a) = In x (1c) 


For sufficiently small values of w, the system's amplitude and frequency characteristics can be represented 
by series expansions: 


2 (28) 

1 k 
= (O)eo"- 2b 
2 (2) 


We investigate the structure of the component of the formula, K (jw) = M(w)eJ?(), in the neighborhood 
of the point w = 0. 


It follows from the analyticity of the system's transfer function in the neighborhood of the point p = 0 that 


K (jo) = 


a* 


In accordance with Formula (1a), this allows us to write the following equality: 


k =0 k=-0 k=-0 


where 


. 


In order to express the coefficients o},, by means of the coefficients ¢, = * @*) (0), we apply Leibnitz' 
formula to the right side of the equation: 


at 


This allows one to write that, for k > 0, 


Turning back to Expression (3), we notice that 


K)(0) = + (5) 


- Analysis of the right side of this equation shows that it is the expansion of the determinant 


My (—1) My Me (—1)k M, 


by the elements of the first row. 
2. Computation of the coefficients o), and G;, is essentially simplified if it is borne in mind that, for odd 
values of k, My = + m‘*) (0) = 0 and, for even values of k, including k = 0, 9 = = A*)(0) = 0. 


SUMMARY 


Expressions were obtained relating the error coefficients with the Cauchy coefficients for the system's am- 
plitude and phase characteristics. 
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REVIEW 


ELECTRICAL DEVICES FOR SOLVING ALGEBRAIC EQUATIONS 


N. N. Mikhailov 


(Moscow) 


The existing root finder circuits are reviewed, classified 
and evaluated from the point of view of the possible automati- 
zation of the root finding process. 


INTRODUCTION 


Many approximate methods for determining numerical values of roots of higher order algebraic equations 
are extremely labor consuming, hence, for the purpose of speeding up and facilitating the finding of roots, com- 
puters which can be called root finders are being developed in many countries. 


The most perfect of these devices find roots automatically and display the results in the form of a root dis- 
tribution picture on a complex plane. The screen of a cathode-ray tube usually serves as the complex plane in 
such root finders; the roots appear on the screen as luminous points, whose coordinates correspond to those of the 
roots, 


’ By observing the root distribution picture of a characteristic ACS (Automatic Control System) equation, 
one can select coefficients which would ensure stability and certain qualitative indexes, such as the degree of 
stability, tendency to oscillate, etc. It is also relatively simple to find by means of root finders the boundaries 
of stability both with respect to one coefficient and to the plane of two coefficients. It is possible to observe 
the migration of roots and traces their trajectory by changing the coefficients of the equation set up on the root 
finder, 


Besides their basic purpose root finders can be used for evaluating polynomials with scalar and complex 
variables, for the analysis and synthesis of Fourier series, for the solution of transcendental equations whose terms 
are represented by power series, etc. It is easy to make up on the basis of root finders devices for conformal 
mapping, used in designing automatic control systems. 

There exist mechanical, hydraulic, electrical, electrolytic and other root finders. Only electrical devices* 
are dealt with in this article, since they are at present the most perfectly developed. 


I. General Conditions 


In a general case roots of an algebraic equation 
w(p)=4,p" +a, a,p’+...+ ap +a =0 (1) 


can be represented in a polar or cartesian system of coordinates in one of the three forms of notation: 


the exponential form 
(2) 


* As far as the author knows one of the first electrical root finder circuits was developed by K. V. Filatov [1]. 
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the trigonometrical form 


P = e (cos 6 + ; sin 9); (3) 


the algebraic form 


p=z+yy. (4) 
When above expressions are introduced into the original equation three forms are obtained for the poly- 
nomial w: ; 
w= +. +a 4... + + ag, 
n n (5) 
w= aye” cos v0 + >} sin = U + jV, 
v=0 


v=0 


(6) 
w= ay +... +4, +... +41 (@+ jy) + do = U (a, y) + jV (wy). (7) 


Let us call Form (5) exponential, Form (6) trigonometrical and Form (7) algebraic, The coefficients as 
well as the roots of the polynomial can be complex quantities, 


The basic principle of operation of electric root finders consists in selecting such values for p and @ (or 
x and y) which convert polynomial w into zero. In the trignometrical and algebraic forms of expression both 
the real part U and the imaginary part V of the polynomial w must become zero simultaneously, Root selec- 
tion can be carried out either by the method of successive approximations or by inspecting the plane of the roots, 


The method of successive approximations consists in changing alternately coordinates p and 6 (x and y) 
in such a way as to decrease the modulus of polynomial w with each change. The coordinates corresponding to 
the zero value of w are those of the root. 


The evaluation of the root by successive approximations can be carried out either manually or automatically. 
In the manual method the root finder operator, while observing quantity | w|, changes the coordinates by hand 
striving to reduce |w] to zero. In the automatic method the operator does not take part in the evaluation pro- 
cessand the required changes of the coordinates are obtained by means of feedback between w and the coordi- 
nates p and @ (x and y). In this case the root finder becomes a closed automatically controlled system. 


The method of root plane inspection consists in moving variable p over the complex root plane along a 


fixed trajectory covering all the points of the area under investigation, and in clamping those values of p which 
convert w into zero. 


The evaluation of roots by this method can also be carried out manually or automatically, In the manual 
method variable p usually describes on the root plane (p-plane) a circle, with its center at the origin of coordinates, 
the radius of which the operator has to change by hand, To this circular movement of point p corresponds the 
movement of point w(p) over plane w along a closed trajectory known as the polar curve. When p passes through 
a point which corresponds to a root of a given polynomial, the trajectory w(p) crosses the origin of coordinates 
on plane w. In the manual method the operator must, therefore, observe the polar curve and changeradius p so 
as to make the curve pass through the origin of coordinates, The corresponding value of the radius is modulus p 
of the required root. Root finders in which the evaluation of the roots is done manually by the inspection method 
can be called semi-automatic. * 


In the automatic method the point determining variable p is given a continuous spiral motion around the 
origin of coordinates in plane p without any assistance from the operator. Usually a part of the spiral trajectory 
limited by a minimum and maximum radii is chosen so that the investigated area of plane p has the shape of a 
ring. Although the spiral trajectory does not pass through all the points of the ring area, it is possible to increase 


* In foreign literature such root finders are often called isographs. 
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the density of spiral winding to such an extent that, with some of the spaces remaining insensitive to the zero in- 
dicator, all the roots lying in that area will be clamped. 


In electrical root finders of any type, owing to their physical limitations, unity is usually taken as the upper 
limit for the modulus of variable p. The lower limit is determined by accuracy requirements and {s often taken 
as 0.1, For determining roots lying outside above limits the equation is transformed by changing variable p. The 
following substitutions are used: 


P (8) 


where k is the constant coefficient and S the new variable. 


In order not to complicate the root finder circuit, the coefficients of the equation to be solved must not 
exceed unity, whichis easily attainable if the original equation is divided by its largest coefficient. 


Fig. 1. Selsyn root finder circuit. 


Some of the root finders are designed to solve Eq. (1) multiplied by p, i. e., without a free term. This is 
done with the object of unifying circuit components. The multiplication by p does not change, of course, the 
roots of the equation, but simply adds one zero root which can be subsequently omitted. In describing concrete 
root finders below it will not be stated whether they solve equations with or without a free term. 


An analytical examination shows that the method of inspecting a complex root plane has a number of ad- 
vantages over the successive approximations method. Thus in manually operated root finders the evaluation of 
roots by the inspection method is much faster than by the approximation method. The actual working of an 
automatic root finder based on the inspection method at the Institute of Automation and Remote Control of the 
Academy of Science, USSR has shown that it can deal with a wider range of problems including the finding of 
roots of algebraic equations. 


In electrical root finders the terms of polynomial w are usually represented by tensions. The characteristics 
of the calculating part of the instrument which produces s these tensions are determined by the form in which the 
polynomial w is presented, whether it is (5), (6) or (7). The basic classification of root finders, therefore, will 
depend not only on the methods and means of finding roots, but also on the polynomial form which the calcula- 
ting part of the instrument handles. By the latter characteristic the existing devices can be divided into three 
types: 1) exponential form polynomial root finders, 2) trigonometrical form polynomial root finders and 3) alge- 
braic form polynomial root finders.* 


* Algebraic form polynomial root finders are usually employed as analog computers adapted for solving equa- 
tions [23, 24]. Hence, they have been excluded from further consideration. 
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Il. Exponential Form Polynomial Root Finders 


In the early electrical root finders, based on the exponential form polynomial w, the evaluation of roots 
was carried out manually by the successive approximation method. Similar root finder circuits were proposed 
by Hart [2], Schooley [3], Elizarenkov [4] and Ginsburg [5]. As an example let us examine Elizarenkov's circuit 
(Fig. 1). 


The selsyn rotors Ly, Lg, Ls, . . ., Ly are interconnected by means of a gear assembly in such a manner that 
when rotor Ly is turned over an angle 6, rotor L, turns over 20, etc, Rotor L, turns over an angle of n@. The 


selsyn stators are fed from a three-phase supply of frequency f. The rotor windings are connected to the modulus 
setting potentiometers P;, Pa,..., Pp. 


The relation between the output tension and position of the slider on potentiometer P, is linear, on P, quad- 
ratic, on P,, of the nth power. The potentiometer slides are connected to a common axis turned by means of 
handle p over equal angles, The slides are connected to linear coefficient setting potentiometers ay, ag, . . ., an. 
The circuit also contains a transformer whose secondary winding Ly is connected to the a, coefficient setting poten-~ 
tiometer, The circuit is pre-set in such a way that, in the zero position of handle @, all tensions in the selsyn 
rotors and transformer winding Ly are in phase and equal in amplitude. 


Tensions at the output of the potentiometers ao, a3, . . ., 4p, in any arbitrary position of handles @ and p, 
can be expressed by means of the symbolic method, well known in electrical engineering as; 


The sum of these tensions Ey can be expressed as 


By + aype?® + aq) 
(10) 


Tension E,, equals zero when +... + +a9=0 4. e,, when p and assume values of 
Eq. (1) root coordinates. The summation of tensions E,, Ey, . . ., Eq is carried out by connecting the coefficient 
potentiometers in series. The zero value of the sum is observed by the operator on voltmeter V. Tumbler switches 
Ty to Tp serve to place signs before coefficients a, to an. 


The differences between this circuit and those mentioned in [2, 3 and 5] are unessential. Instead of selsyns 
Hart and Schooley used ac generators. In Ginsburg's instrument the polynomial w is set up in the calculating part 
of the device in its exponential! form, but is later divided into its real and imaginary components, which provides 
an opportunity of representing w as a point on a cathode-ray tube screen, simulating plane w. 


The accuracy of root evaluation in the selsyn circuit is 2-10% for the modulus and 3-5° for the angle, The 
sources of error include, in addition to the variation in circuit components, distortion of the carrier waveform and 
unbalance between the selsyn phase supply voltages. 


In 1955 Parker and Williams [6] constructed on the basis of the same circuit a more accurate instrument, 
using high quality components (magslips, spiral potentiometers, summation amplifiers) and a special circuit for 
eliminating supply phase unbalance and carrier distortion, which reduced the error to fractions of one percent. 
The gear assembly was considerably simplified by making some of the selsyns rotate in the opposite direction 


to the others. The powers of p were obtained by a cascade connection of linear potentiometers, separated by 
amplifiers. 


Bubb [7] developed a root finder circuit essentially different from the above, His circuit is based on linear 
potentiometers and operational amplifiers, i. «, high gain amplifiers with a large negative feedback. As is known 
the input voltage Ej and output voltage E, of such amplifiers bear the relation 


(11) 


where z, is the input impedance and zg the feedback impedance. Depending on the nature of z; and zz Expres~ 
sion (11) becomes; 


Fig. 2, Operational amplifier root finder circuit. 


a) with z; = zg =R 
(12) 
4 4 
with 2; R, % 
Ey = (14) 


In addition to above operations the amplifier can also register a sum. Utilizing these properties of opera- 
tional amplifiers Bubb designed a simple circuit which multiplied ac tensions by complex quantities p = x + jy. 
It is easy to show that tension E, at the output of element 1 (Fig. 2) is relatedto tension /,=|£,\e?"/' at the 


input by the expression 
BE, = (x + jy) Ey = pEo. (15) 


The values of |x| and|y| are set on linear potentiometers Px and Py and their signs are determined by 
the position of tumbler switches T, and Ty. 


r 
r 
d 
F 
s 
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Fig. 3. Operational amplifier and magslip root finder circuit. 
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If several such circuit elements are connected in series then the tension at the output of the first element 
will be Ey = pEp, of the second —E, = p’E,, the third Es = p’E9, etc. All the slides P, are connected to a common 
axis |x| , and the slides Py to axis] y| . The corresponding tumbler switches must "alto be mechanically inter- 
connected. A similar circuit element is used for multiplying the obtained powers of p” by complex coefficients 
a, =b, + jcy. Next the terms a, p” Ey are added up on the operation amplifier Y5, whose output tension is mea~ 
sured in instrument V. The roots are evaluated manually by the successive approximations method. This circuit 
formed the basis of Lukashevich's root finder [8]. 


A disadvantage of setting up variable p in Cartesian coordinates is the switching required to pass from one 
quadrant to another. This disadvantage is eliminated if polar coordinates are used for the solution of the problem. 


Such a changeover is easily carried out [9] by replacing the potentiometers P, and Py by magslips M as it is shown 
in Fig. 3. 


The magslip cosine winding is connected to the operational amplifier through a resistance, whereas the 
sine winding is connected through a capacity; thus the circuit consisting of a modulus p setting potentiometer, 
a magslip and one operational amplifier provide for the multiplication of the input tension by p = p(cos 0 + 
+ jsin@), in which 6 can vary between 0 and 360°. The switch T serves to reverse the sign of the tension. 


III, Trigonometrical Form Polynomial Root Finder 


Basic operations which must be performed by a trigonometrical form polynomial root finder when setting 
up terms a, p’ cos v@ and a,p” sinv@ of Expression (6) are; obtaining powers of p, taking into account trigo- 
nometrical relations of the type of cos v@ and sinv@ and multiplication. For multiplication electromechanical 
devices are mainly used (potentiometers, variacs, magslips, etc.) in which one of the factors is the electrical 
tension and the other mechanical displacement. Let us call multiplicand the factor provided electrically and 
multiplier that provided mechanically. Depending on which function serves as the multiplier the trigonometrical 
form polynomial root finders can be divided into three groups: in group I p” serves as the multiplier, in group II 
it is cos y@ and group III includes root finders with a purely electronic method of multiplication. 


Fig. 4. Sine-cosine potentiometer root finders, 


Group I circuit. The first circuit of a trigonometrical form polynomial root finder was proposed by Eliza- 
renkov [4] (Fig. 4). It consists of sine-cosine potentiometers P, to Py, whose slides are connected by gear assem-~- 
blies in such a way that by turning with handle @ the first potentiometer over an angle 0, the second is made to 
rotate over an angle 20, the third over 30, etc, The inputs of the potentiometers are connected to a common 
de tension Ey. The sine and cosine outputs of potentiometer P; are connected to two similar modulus setting 
potentiometers with a linear winding, P, is connected to potentiometers with a quadratic winding, Ps to potentio- 
meters with a cubic winding, etc. The slides of all these potentiometers are connected to a common axis for 
ey p, and the potentiometers themselves to the linear coefficient cere potentiometers aj and af, ag and 


aj, aj and aj, etc. At the series connected outputs of aj, aj, ag, . . . and af, af, af, . . . two respective voltages 
n n 

are formed U = >! ap” cos v0 The evaluation of 
vel vel 


roots is made manually by the successive approximations method. 


D 


Root finding in this circuit was made automatic by G, M, Zhdanov [10] who proposed to operate the p -axis 
by means of a servo-unit whose input is connected to tension U and to connect tension V to the input of a servo- 
unit controlling the @-axis. The servo-units stop when tensions U and V become zero, i. e., when p and @ be- 
come coordinates of the required root. As it has already been pointed out such a root finder becomes a closed 
automatically controlled system, It would appear that the stability of such a system and the iime taken to achieve 
balance will depend on the character of coefficients in the original equation, No experimental data for the opera- 
tion of such a root finder exists, since as far as it is known the instrument was never constructed. 


If axis 6 in the circuit of Fig. 4 is rotated at a constant speed roots can be found by the inspection method, 
This rotation will also make harmonic frequencies appear at the outputs of the sine-cosine potentiometers. It ts, 
obviously, better to generate such frequencies in a purely electronic manner by means of a harmonic oscillator, 
which was in fact carried out by Tischner [22] and Bader [11] (Fig. 5). The formation of harmonics in this case 
is carried out by means of generator G, a frequency multiplying circuit FMC and filters F; to Fp. 


The use of two modulus and two coefficient potentiometers in each channel is a disadvantage of the sine- 

cosine potentiometer root finders, In the harmonic generator circuit each channel has only one such potentio- 
meter and the formation of the sine and cosine components of Expression (6) is achieved by phase splitting ele- 
ments PS. 


FMC 


Fig. 5. Harmonic generator root finder. 


The aggregate tensions U and V appear on the horizontal and vertical deflecting plates of a cathode-ray 
tube CRT, whose screen serves as plane w on which appears the polar curve. The operator turns handle p striving 
to make the curve pass through the origin of the coordinates. The value of p corresponding to this position is 
the modulus of the required root. For determining angle @ on the polar curve a mark {s formed electronically 
(in the shape of a bright spot or a break in the curve) which represents a point on the circumference of radius p 
on the plane w. The operator controls the position of the mark and observes the value of 6 on a scale of the 
measuring circuit MC. Having placed the mark on the polar curve point lying at the origin of the coordinates 
on plane w (the screen of the tube), the operator reads off the scale the value of the angle of the root. 


Choudhury [12] proposed a root finder circuit using an artificial long line. As is known tension E at a 
distance x away from the end of such a line is equal to 


= Bch yy + 1,2 sh (16) 


where &, and }, are the tension and current, respectively, at the end of the line, z is the characteristic im- 


pedance and y = 6 + ja is the propagation constant. In a line without losses 8 = 0 and Eq, (16) becomes 
= cos ay + sin ay. (17) 
In the case of free running with open circuited output terminals we have 


E = i, cos ay. (18) 
* ‘cosh’, ‘sh’ = ‘sinh’ — Publisher. 
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With short circuited terminals we obtain 


=jl,z sin ay. (19) 
If a long line consists of nJ similar low pass filter sections, x denotes the number of sections between the 


end of the line and a given point. Fora point v! sections away from the end of the line expressions (18) and 
(19) will become 


£,=£,cosval, it, = jl, sin val (20) 
or if one assumes al = @: 


cos v6, = jl,z sin we. (21) 
In the passband of low-pass constant-k type filters * a= 2 arcain £ , where f, is the cutoff frequency. 
c 
Hence, 


= 2laresin (22) 


In Choudhury's root finder (Fig. 6) the artificial delay line DL is connected to the swinging frequency gen- 
erator SFG, controlled by the sawtooth oscilloscope time-base generator SG, The limits of frequency f varia- 
tions are selected in such a way that 0.383 = f/ f, =< 0.707, and 7 is taken as 8, © varies between 0 and 360°, 


Fig. 6. Delay line root finder circuit. 


Tension E,, is taken from the output of the required filter section, connected through a separating amplifier 
to the coefficient setting potentiometer a,,and then multiplied byp” and simultaneously added toother tensions. 
In this circuit, as distinct from the circuits described above, it was possible to use linear potentiometers for multi- 
plying by various powers of p. The output tension E,, corresponds either to component V or U of Expression (6) 
if the output terminals of the line are closed oropen, This tension is impressed on the vertical plates of the os- 
cilloscope whose horizontal time-base controls the frequency f. 


* Filters with their series and shunt arms consisting of two-terminal networks whose product is a constant are 


known as “constant-k" type filters. 
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A rough evaluation of roots is supposed to be done manually by the inspection method alternately shorting 
and opening the line output terminals* with the switch T, An accurate value of coordinates is obtained by the 
successive approximations method. 


Group Ii circuits. The magslip circuit [14-16] is being widely used (Fig. 7). Magslips MG permit one to 
vary rapidly and thus change over from the successive approximations method [14] to the manually operated in- 
spection method (15, 16]. Continuous variation of coordinate 6 is achieved, in this case, by means of a motor 
M which turns the magslip rotors through a gear assembly. The use of magslips made it possible to change the 
circuit over to ac working. 


Fig. 7. Magsiip root finder circuit. 


Various powers of p are obtained by means of a cascade connection of linear potentiometers with a com- 
mon slider axis. In order to prevent the possible loading of potentiometers, they are separated by amplifier with 
a gain equal to one, The values of coefficients are set on linear potentiometers a,,. The formation of the real 
U and imaginary V components of Expression (6) is obtained from the carrier tension by grouping in a series con- 
nection the cosine and the sine windings of the magslips M. The evaluation of roots is carried out in the same 
manner as in the harmonic generator circuit. 


Adler's instrument [17] is a modification of above circuit in which the magslips and modulus setting poten- 
tiometers are replaced by tapped transformers; transformers are also used for setting up coefficients. 


Calvert,Johnston and Singer [18] proposed an automatic magslip circuit in which a motor drives with con- 
stant speed not only the @ axis but also the p axis. Variable p = pe j9 describes a spiral p = (Wp /we)8. The 
working accuracy depends on the density of the spiral winding, hence the ratio wp /ug was taken me 0.01, and 
the frequency f as 400 cps. The period of p variation in the interval between 1 and 0.1, 1. e., the inspection 
time of plane p ring area, bounded by radius pmax = 1 and pmin = 0.1 is equal to 18 minutes. 


* The method of alternately shorting and opening the end of the line cannot be considered convenient. In order 
to generate U and V simultaneously two similar delay lines are required, which is not expedient. It is possible, 
however, to indicate a delay line root finder design free from above disadvantage. 

If the line is terminated in its characteristic impedance and it is noted that I = E,,/z it becomes possible 
possible to write Expression (17) in the form 


E, 
=k, cos ay +- by — kt, ei” (23) 


If E,, is multiplied by a, and p” and added to the remaining tensions, the tension E, obtained at the out- 
put of the circuit will represent the entire polynomial w in the form (5) and not a part of it only as is the case 
in Choudhury’s circuit. The circuit thus obtained should, of course, be classified as an exponential form poly- 
nomial root finder. 
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The demodulated tensions U and V are impressed on an electronic zero indicator which clamps the zero 
values of | U| and |V|,stops the motor and informs the operator by a signal. 


An accurate evaluation of the root thus obtained is made manually by the method of successive approxi- 
mations, In this root finder certain measures have been taken to improve accuracy. For instance, the gain of 


the separating amplifiers changes automatically according to the value of p, helical potentiometers are used 
for setting moduli. 


The main object of making root finding automatic is not so much for the purpose of saving time in solving 
equations as for obtaining (usually on a cathode-ray tube screen) a picture of root distribution which would faci- 
litate the observation of root migration when any coefficient in the equation is changed, For this purpose the 
inspection process of the complex plane selected area must be repeated periodically, and the frequency of repe- 


tition Q must be sufficiently high in order to make it possible to obtain root trajectories in the form of curves 
on the tube screen, 


The above magslip automatic circuit does not permit one to obtain root trajectories owing to the very low 
w 
repetition frequency (Q = rg = 107° cps), moreover, the time spent on root finding is probably increased instead 


of being decreased as compared with non-automatic circuits [15-17]. The root finding can only be speeded up 
in this circuit by increasing the carrier frequency f, limited by the magslip operating supply frequency, which 
at the present time cannot exceed 1000 cps. Another drawback of all carrier automatic root finder circuits is 

the necessity of demodulating tensions U and V before they can be impressed on the zero indicator. In order not 
to distort the envelope in demodulation, the repetition frequency Q has to be kept low. The use of a carrier fre- 
quency, therefore, in automatic root finder circuits must be considered inexpedient. 


A root finder without a carrier on the basis of a manually operated inspection method was constructed by 
Marshall [19] (Fig. 8). The powers of p are formed in it in the same manner as in the magslip circuit, the two 
differ, however, in the fact that the former has its potentiometer cascade fed by dc. Amplitude modulated square- 


w 
wave pulses of a frequency v = are obtained from tension Ego” by means of an electromechanical modulator . 


The modulator contains n rings mounted on a common axis; the rings have equally spaced conducting and insu- 
lated segments. The first ring has one conducting segment, the second has two, the third three, etc, The conduc- 
ting segments are connected to a common point of the circuit, so that when the brush is in contact with one of 
these segments the tension at the output of the modulator is zero, when the brush comes over an insulated seg- 
a ment the tensions becomes equal to Eygp”. When the 
motor M rotates at a constant speed the y -th channel 
modulator will produce a square-wave tension whose 
amplitude will be proportional to p”, Filter F, picks 
out of these oscillations the first harmonic tension which 
is multiplied on a potentiometer by coefficient a, and 
then impressed on a phase-splitting RC network whose 
output tensions are equal to Egay p” cosyvWg t and Epa, - 
p” sin vwet. 


The signs of the coefficients are set by the switch 
T. The aggregate tensions U and V form on the screen 
of the cathode-ray tube a polar curve. The quantity 


1s mace to equal 60 cps. 


Lofgren [20] substituted the electromechanical 


Meas, 


L 


eitguit | modulator by an electronic one, this allowed him to 

rit pss raise we and change over to an automatic method of 
root finding, Automatic variations of p can be produced 
without any moving parts by utilizing the exponential 
relationship of p to time, Such a law of variation of p 
can be relatively easily produced by periodic switching 


Fig. 8. Square-wave modulator root finder circuit. 


of RC circuits. As is known the raising to power vof an exponential quantity will produce another such quantity 
but with a time constant 1/y times the former value. If the time constant of the RC networks r = RC ts selected 
in such a way that among the various channel time constants the following relation holds: 


the tensions produced by the networks of the second, third, etc. channel will be respectively the square, cube, etc. 
of the first channel RC network tension. The switching frequency in Lofgren's root finder is taken as Q = w)/2m = 
= 1 cps and the value of 1, is selected to make p vary in the limits 4/, <p <1. 


The square-wave harmonic frequency oscillations are generated by a special frequency dividing circuit. 
The quantity ra is approximately equal to 1000 cps. The formation of the sine and cosine components of Ex- 


pression (6) is achieved without phase~splitting networks by an operational amplifier-summator and summator- 
integrator. 


The picture of root distribution of the required equation is produced in the form of light spots on a cathode- 
ray tube screen with a long afterglow. The screen should simulate part of the complex root screen. For this pur- 
pose the tube deflecting systems are supplied with tensions which displace the ray along the same spiral p = 
= 9", as the one controlling the automatic variation of p in the counting part of the instrument, Normally 
the ray is suppressed and it is only switched on for the instant when the tensions U and V, which represent the real 
and imaginary components of the polynomial, simultaneously become zero when they are clamped by the zero 
indicator. 
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Fig. 9. Electronic multiplier root finder circuit. 


For an accurate evaluation of the root coordinates thus obtained a special measuring circuit is used. This 
circuit produces on the tube screen a movable mark, whose displacement can be controlled by handles with scales. 
When the mark is made to coincide with a root which appeared on the screen its coordinates can be read off the 
scales, 


The automatic root finder of the Automation and Remote Control Institute of the Academy of Sciences, 
USSR was constructed on the basis of the square-wave modulator circuit [21]. 


Group III circuits. The first automatic root finder circuit which displayed the root distribution on a screen 
of a cathode-ray tube was developed by Rasch [22] who proposed a purely electronic method for forming p” cos v- 
‘Wet and p” sinvwet components, This method is based on the fact that when oscillations of a fundamental 
frequency pcoswet are impressed on a nonlinear circuit, at the output of that circuit among other components 
will also appear tensions p” cosywet which can be filtered out. 


Rasch's circuit (Fig. 9) consists of a generator G which supplies the fundamental frequency, adjusting unit 
A, for adjusting the amplitude of p, and the above-mentioned circuit with a nonlinear characteristic NLC which 
serves simultaneously for multiplying and raising to the power. This circuit is followed by filters Fy to Fy which 
pick out the required frequencies and by units for setting coefficients a,,. The separation of tensions into their 
sine and cosine components is achieved by means of phase splitters PS. 


The picture of root distribution is obtained in the same way as in the automation square-wave modulated 
root finder described above. 


On the basis of this circuit Glubrecht [22] constructed an automatic inspection method root finder with an 
alternative manual operation; for the latter purpose a screen representing plane w was incorporated in the instru- 
ment in addition to the screen simulating plane p. The change over to manual operation is intended for accurate 
evaluation of root coordinates. For generating oscillations of the form p” cos vw gt instead of a nonlinear circuit 
an electron tube is provided whose control grids have tensions p cos wet and p”~! cos(y — 1) wp t impressed on 
them. The output tension contains component p” cosvw, t which is filtered out. 


IV. Root Finders for Solving ACS Characteristic Equations 
An automatic control system (ACS) characteristic equation can be given in the form 


G(p)+1=0. (25) 
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Fig. 10. Circuit for solving ACS characteristic equations. 


If the transfer function of an open system G(p) is given by the expression 


+ _ P™ + by 


G(p)=k 
a, + +... + 4% (26) 


after substituting G(p) in (25) by (26) and eliminating the common denominator we obtain 


+. +b, + bo) = 0. (27) 


By grouping similar terms the equation is reduced to the usual (1) form of an algebraic equation which can 
be solved on the types of root finders described above. If, however, the object is to investigate the relation between 
the initial coefficients a,,, b; and k and the position of roots, the reduction of the equation to Form (1) fs not 
expedient, since this only complicates investigation. It is desirable to be able to solve the characteristic equa- 
tion in the form (27). 


+ Ve, 


It ts easy to construct devices for the solution of such equations on the basis of root finders of any type. For 
this purpose it is necessary to connect in parallel with the a, potentiometers others for setting the by coefficients 
and to connect them to a circuit similar to the one to which the ay potentiometers are connected. 


As an example of such a device a circuit based on a selsyn root finder is given in Fig. 10 (secondary details 
have been omitted). Tensions E, and £} at the outputs of the summing amplifiers can be written down as 


hed m m 
ved veo 


Fig. 11. Circuit for multiplying tensions by (p—q) and (p* + 3p + ») 


Coefficient k is set on the slider of potentiometer P,. The output tension measured on meter V is equal to 


> ap’ +k (29) 


which corresponds to Expression (27). ; 
The transfer function G(p) is more often given in the form 


(p — 21) (p — 22)... (P—2,,) 
(p— 41) (P—q2)---(P— 


G (p)=k (30) 


where 24, Za, . . ., Zp are zeros and qy, Gg, . . ., Gn are poles. Usually one is interested in the relations between 
the zeros and poles, and the position of roots of the characteristic equation, therefore, the transformation of 
Expression (30) to the form (26) by the opening of brackets is not admissible. 


The construction of a root finder for the solution 
of ASC characteristic equations when G(p) is given in 
the form (30) is a more complicated problem mainly 
because products of complex numbers have to be formed 
in the counting part of the instrument. Nevertheless, 
by taking as the basis a root finder with operational am-~- 
plifiers and magslips it is possible to make up a circuit 


which would evaluate roots of such equations. 


Fig. 12. Block schematic of a root finder for the solu- 
tion of equations : The substitution of G(p) in (25) by (30) and the 


elimination of the common denominator will give 


(p= 1) (P= 92) - (P — 21) (p —22)--- 


Tensions E, and Ep: 
= Eo(p —%) (Pp — 92) ---(P~ = 
= Ko(p— 2) (p — 2x) -- -(P— 2m). 


k(p—2)(p— 22) .. —2,,) = 0. (31) 
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The poles q,, and zeros z; can be both real and complex conjugate numbers, Hence, the circuit must con- 
tain components capable of multiplying tensions by complex expressions of the type of (p— q) and (p* + 5p +n) 
where q, 6 and p are realnumbers. For these operations circuits shown in Fig. 11 can be used. The output ten- 
sion in the circuit of Fig. 11a is 


Eo-= (p— (32) 


The circuit in Fig. 11b reproduces the relationship 


(33) 


The signs of coefficients qg, 6 and » are produced by switch T. 


A circuit for producing the left-hand side of Eq. (31) and consisting of such components fs shown in Fig. 12. 
Amplification factor k is set by means of potentiometer P,. A root finder constructed on the basis of this circuit 
can be made to evaluate roots either manually or automatically. 


SUMMARY 


1. It is advisable, according to the nature of the problems to be solved, to have two types of root finders; 
semi-automatic and automatic, The former are intended to find roots and the latter mainly to investigate the 
relation between the roots and the coefficients of the equations, 


2. It is advisable to construct automatic root finders in the form of open systems which find roots by the 
inspection method. 


3. The use of a carrier in root finder circuits limits the repetition frequency and thus makes it difficult 
to obtain trajectories of root displacements, 


4. It is possible to form a polynomial in a trigonometrical form without a carrier, hence, the use of this 
polynomial form permits one to construct the fastest automatic root finders. 


5. The limitations placed on the carrier and hence on the repetition frequency by the use of selsyns and 
magslips are absent in the delay line circuit. Hence, if a polynomial is formed in an exponential form (and not 
a trigonometrical form as suggested by Choudhury) this circuit can serve as a basis for the construction of a fast 
root finder, 


6. The existing types of root finders are not suitable for solving ACS characteristic equations, it is, there- 
fore, desirable to construct devices which would investigate the relations not only between the roots of a closed 
ACS equation, and the poles, but also the zeros of an open system transfer equation, Possible circuits for such 
devices have been suggested in part IV of the article. 


1. By utilizing the principle of continuous inspection (scanning) of a complex plane it is possible to con- 
struct, on the basis of existing root finders, a device for conformal mapping. 
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CHRONICLE 


NATIONAL CONFERENCE AT BUCHAREST ON THE QUESTION OF 
AUTOMATION OF PRODUCTION 


During June 5-8, 1957, the Second National Conference on the Question of Automation of Production was 
held in Bucharest (Rumanian National Republic). Invited to participate in the Conference were scientists and 
specialists in automation from Bulgaria, Poland, Hungary, Rumania, Czechoslovakia, the USSR and Yugoslavia. 


In the majority of cases, these were contributors to academic, research and professional institutes who 
simultaneously engaged in pedagogical work. Among them were: Corr. Mem, Bulgarian Acad, Sci. A, Bolev- 
skii, Prof. D. Mitrovich (Yugoslavia), Dr, I. Benesh (Czechoslovakia), Prof. Libson (Poland) and others, 


A total of 99 papers were presented at the Conference, including 91 from Rumania, 2 from Czechoslovakia, 
2 from Hungary, 1 from Poland, 2 from the USSR and 1 from Yugoslavia. 


The delegation from the USSR included two contributors to the IAT AN SSSR, Prof. N. N. Shumilovskii and 
Eng. I. M. Makarov. 


Besides participating in the work of the Conference, our delegates had the tasks of familiarizing themselves 
with the arrangements for research work in automation and related fields in the scientific institutions of the AN RNR 
and in the professional institutes, familiarizing themselves with the state of automatic production in the factories 
and, in addition, establishing scientific contacts with the scientific institutions and the scientists of the RNR and 
of the other countries represented at the Conference. 


The National Conference on Automation was called at the initiative of the Academy of Science of the RNR. 
The Academy does not have an institute devoted to the questions of automation, but there is a committee on 
automation, headed by Academician I. S, Georgiu, The Committee has two scientist-secretaries, contributors 
to academic and professional institutes, to institutions of higher learning, industrial plants, ministries and Gosplan 
RNR. All the organizational work involved in summoning the Conference was done under the aegis of this Com~- 
mittee, 


Several months before the opening of the Conference, the organizing committee sent requests to the de- 
sired foreign participants requesting that texts of papers be furnished by the proper dates. 


The Conference opened on June 5 in one of the buildings of the Bucharest State University. Presiding Aca~ 
demician I. S. Georgiu opened the session with a few words welcoming the guests and participants in the Confer- 
ence and outlining the problem which the Conference was to consider, Further words were added by Prof, N. N. 
Shumilovskii, welcoming the participants in the name of the Institute for Automation and Remote Control of the 
Academy of Sciences of the USSR, and wishing success to the Conference in its work, Short welcoming addresses 
were made by representatives of the other countries, after which the order of business of the Conference was an~ 
nounced, 


The participants were given programs containing the theses of the Rumanian delegates" papers, these pro- 
grams being printed in the Rumanian language, which somewhat attenuated the intensity of the discussions. At 
the plenary session, papers were presented by guest delegates Prof. N. N. Shumilovskii, Dr. 1. Benesh and Prof, 
D. Mitrovich, dealing with the state of automation and measurement technology in their countries in various 
domains of industry. Then followed Rumanian delegates; Academician G. K. Moisil, Prof. K. Penesku, Engr. 
D. Damsker, Engr. M. Maresh, Corr. Mem, RAN A. Avramensku and others, who, in their papers, dwelt on the 
questions of automating various branches of industry in the RNR, on questions of terminology in automata, and 
also on questions of preparing personnel for the automation of industry, 


Further work of the Conference was carried out in three sections. The first section dealt with the theoretical 
basis of automation. The second section considered the general questions of automation, measurement and con- 
trol of productive processes, and the third section was concerned with questions of automation, measurement and 
technological control in various branches of industry. 


The first section considered questions of linear and nonlinear theories of automatic control, the theory of 


relay~contact designs, statistical methods of investigating automatic control systems, and certain questions in the 
theory of electric drives. 


Among the papers presented at this section, we mention the following: “Synchronous Motors and Motors 
with Variable Inductance as Servomechanisms with Proportional Control” by Corr. Mem, AN RNR M. Marinesku 
and Engr. G. lankylesku, "A Contribution to the Simplification of the Stability Conditions of A. I. Lur'e* by 
Engr. V. Popov, “Construction of the Transient Response in Systems with Nonlinear Regulators" by Engr. D. Dam- 
sker, “Active Correcting Links" by Engr. K. Vazak, “A Method of Simplifying the Analysis of Continuous Linear 
Systems of Automatic Control” by K, Penesku, "A Method of Designing Discontinuous Automatic Control Systems * 
by Engr. N. Shtefanesku and “The Characteristic Equation of a Tripping Relay” by Academician G. Moisil of the 
AN RNR, 


Papers presented in the first section by foreign delegates included: “Some Questions in the Theory of Non- 
linear Automatic Control Systems" by Dr. S. Vengzhin (Poland), “Statistical Method for Investigating the Dynamics 
of Control Systems" by Dr. I. Benesh (Czechoslovakia) and "New Method of Determining the Characteristics of 
Single~Phase Electric Spindles" by Engr. L. Ilanoki (Hungary). 


In the papers of the second section, in addition to questions of automation, measurement and control of 
individual productive processes, there were considered the questions of simulating, creation of digital computers 
and also methods for the design and construction of individual elements of automatic devices. Together with 
these, papers were presented dealing with the economic questions of automation. A total of 27 papers were given 
at this section, including: “Economic Considerations in Automation in Capitalist Countries" by Engr. E. Balash, 
“Automation in Countries Within the Socialist Camp” by I. Teodorov and I. Lemin, "The Electronic Computer 
of the Institute for Nuclear Physics of the AN RNR” by Engr. V. Tom, “Electrical Measurement of the Moisture 
Content of Materials" by Engr. M. Steru, “On the Planned Radio-Tele-Control Systems to be Applied to Objects 
Moving at High Speeds” by Engr. R. Konstantinesku and “Questions as to the Method of Designing Determining 
Amplifiers for Electronic Models” by Engr. S. Shekhter. Engr. I. M. Makarov (USSR) gave a paper on the subject, 
"Choosing Optimal Electrodynamic Couplings for Automatic Control Systems.” 


The papers of the third section presented the practical results of the work carried out in the RNR in the 
domain of automation of different branches of industry: petroleum, chemical, light industry, textiles, food, print- 
ing, and others, Some attention wasalso given to work in the automation of communication. That more than 
30 papers on the automation of the national economy of the RNR were given bespeaks the fact that, in Rumania, 
a sufficiently large group of specialists in the domain of automation has been assembled. The plenary session 
of the Conference was held on June 8. Both scientists representing the ministeries and Gosplan of the RNR and 
guests took part in the discussions of this session, In his concluding remarks, Academician I. S. Georgiu men- 
tioned the increased interest in Rumania in the problems of automation and the marked rise in the general level 
of research and planning of automatic devices, which argues well for the further development of automata and 
measuring technology in the Rumanian National Republic. 


In the resolutions of the Conference, together with important organizational and practical measures bear- 
ing directly on the questions of the development of the theory and practice of automation in the RNR, particular 
attention was given to the necessity for increasing scientific and technical cooperation with other countries and, 
as the first step, with the countries in the socialist bloc. 


To this end, it was proposed to develop the interchange of scientific and technical documentation, to work 
out a program of national conferences on question of automation with participation of the other national demo- 
cratic countries, to organize exchange visits of specialists from Rumania and from foreign countries, in particular 
from the national democratic countries, to take steps to include Rumania in the recently created International 
Federation for Automatic Control. 


Not only in the Conference resolutions, but also in personal conversations with the representatives of a num- 
ber of countries, there were expressed hopes for the further strengthening of scientific ties, and of organizing a 
single group of specialists in the domains of automation, remote control and measuring technology. It was pro- 
posed to summon a conference to consider the results of the most important scientific and practical work carried 
out in the countries of the democratic bloc, to interchange regularly printed work so that they might be periodic- 
ally published in these countries, to exchange attempts to apply certain methodologies in scientific investigations, 
and also to organize joint scientific work in the domain of automata. 


The delegates to the Conference had the opportunity to visit the Institute for Nuclear Physics, located not 
far from Bucharest, where a working model of an electronic digital computer was demonstrated, 


On June 9, after the close of the Conference, guests and members of the organizing committee had an ex~- 
cursion to a newly constructed, well-designed petroleum refinery in Ploesti. 


On returning from the excursion, the guests were taken in hand by the scientist-secretaries of the Presidium 
of the AN RNR and held conversations with the leading scientists of the AN RNR. 


In addition to the organizations to which visits were provided for in the Conference program, the Soviet dele- 
gates visited the Energetics Institute of the AN RNR, the Electro-Technology Institute and the Military-Technical 
cs Academy in Bucharest, and also a chemical factory producing superphosphates and sulfuric acid. 


In conclusion, it is desirable to mention that, both on the part of the organizing committee and on the part 
of all its members, the delegation from the USSR experienced the burgeoning of friendly, warm relationships. 


I. M. Makarov and N, N, Shumilovskii 
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